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CHAPTER 1 


INTRCDUCTION 

1*1 Motivation* 

The mathematical modelling of many physically signifi- 
cant problems/ such as phase change in chenical and metalluirgi- 
cal processes, or in glass, plastics and oil industries, or 
in the manufacture of semiconductor materials, leads to boundary 
value problems involving partial differential equations ~ 
especially of parabolic type - with the feature that a part of 
the boundary is not prescribed in advance, but depends on 
certain properties of the solution itself* Usually, the 
d^endent variable r^resents either tenperature, or concentra- 
tion of a solute in a solvent, or pressiire of fluid diffusing 
through a porous medium, or vapour density in an unsaturated 
mixture# These problems are traditionally called free boundary 
problems or more classically Stefan problems# An example of 
such problems is the well-known problem of melting of ice, 
which was first studied by J# Stefan in early 1839, afterwhom 
this class of problems is named* The basic difficulty in such 
problems traces its origin to the determination of a part of 
the boundary in the process of solution of the problem itself* 
This naturally leads to a nonlinear problem, which is always 
harder to solve both analytically and numerically, than the 
corresponding problem with a fixed boundary* 
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In the last few years, several numerical techniques 
like Finite Difference (FD) and Finite Element (FE) Methods 
are being developed for numerical approximatiois to solutions 
of the free boundary problems* However, the finite difference 
methods have been used extensively for numerical solutions of 
these problems (see for recent surveys Purzeland [ 25 ] , 

Nitsche [ 43 ] t [51 ] and Tarzia [59 ])• But in recent years 
several finite element methods have been proposed in view of 
the increased speed and storage capacdty of present day 
Gomputers* As several of the proposed finite element methods 
lack sufficient mathematical justification, we consider it 
worthwhile to present a brief survey of this and make an 
attempt to justify some of them, under suitable regularity 
assunptions* 

Since the dominant difficulty is introduced by the 
unknown free boundary, a basic technique used in gtmeral is 
to have a co-ordinate transformation in order to fix the 
boundary* The penalty one pays for it is the introduction of 
a highly nonlinear term into the transformed problem# Earlier, 
Nitsche using this technique initiated the study of finite 
element error analysis for the single phase linear Stefan 
problems in one space dimension* The present study is a 
contribution in the same direction for the nonlinear case* 

It partially answers the open problem, cited in the workshop 
on Free Boindary Problems held in Berlin 1977 (see Hoffinann 
[29/1], pp* 185)# 
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Essentially we are dealing with two types of problerris^ 
classified according to the form of the transformed system^ In 
the first case, the transformed problem consists of a non- 
linear initial boundary value problem in a fixed domain and 
two initial value problems in ordinary differential equations 
(see Chapter 2# 3 and 4), while in the second case the 
transformed problem is a coupled system of a nonlinear initial 
boundary value problem and two initial value problems in 
ordinary differential equations (see Chapter 5)* For both the 
cases the object is to give an error analysis feor the Galerkin 
methods in different norms • 

1*2 Notations and Definitions * 

In this Section, we introduce some notations and 
definitions which will be used throughout this work* All 
functions treated here are real valued functions* For each 

i 3 f 

integer j >0, D-'f = ■2— i* 

For each t > O, let Q(t) c R be a bounded domain* 

On the space of square integrable functions on Q(t)|. the 
2 

L —inner product and norm are denoted by 

(u,v) = / uv dx 
Q(t) 


and 


llull"^ = (u,u) respectively* 
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As usual, L (Q(fc)) shall denote the apace of bounded measurable 
functions on Q(t) and 

1 lul 1 ^ s ess sup ju(x,t)} < 

L (Q(t)) xPoft) 

The support of a function f : Supp (f) = { x : fix)’ ^ 0}* 

For any nonnegative integer a, cP^(0(t)) denotes the 
space of functions with continuous derivatives upt^o and 
including order a in Q(t), Ag usual, C^(0(t)) denotes all 
c’^(o(‘t)) functions having compact support in oCt), \irhile 
C~(Q(t)) is the space of all infinitely differentiable function 
with compact siopport in Q(t)* 

For each non— negative integer m, the Sobolov space of 
order m on Q(t), denoted by H™(Q(t)) ( respectively, I^(Q(t)), 
is defined as the closure of C^(o(t)) (respectively, C^(0(t)) 
with respect to the norm 

2 ^ i 2 

lluir^ * E / lDJu(x,t)l ^ dx. 

H^(Q(t)) jrO D.(t) 

Equivalently, H'^(Q(t)) is the space of all functions whose 
distributional derivatives of order less than or equal to ra 
are in L^(CJ(t)) - Clearly, L^(Q(t) ) a H°{o(t)) * f^"(Q(t)) is 

a Hilbert space with inner product defined by 

m . . 

(u,v) „ = E / Ir u Irv dx- 

pf'Cmt)) j=OQ(t) 

for a coitplete discussion of Sobolev spaces, see Adams [ l ] • 
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We also define oCt )) for each t > 0 as the 

closure of c'*’(0(t)) with respect to the norm 


m 


i lul ] 


w^'^^CoCt)) j=0 


= r iidJuIi 


L (o(t)) 


In case I =s (0,i) =* 0(t)/ then we shall omit I from 


H^(l), L ( l) and l) and the norm in H^( l) will be denoted 


by I 1 .1 1^. 


If X is a normed linear space with noxm U *1 and 
<P : (a^b) X, then we denote ?ay 


w^'^(a,b;x) 3=0 L^(a|b?x) 


and I I VI I 1 , ^ is accordingly defined* 

w^'”(a,b?x) 


In case (a*b) = (0,T)' and X = h”’ or we write 


' V m I • Vl 1 V rr m ^ V «• m « 


for llvri 




d<f> 


For convenience, we use <P„ *= ^ =* f 

X <jx XX gx*^ 


^t * ^ V(l) = ^(l,t), if 91 = <)>(x/t)* Throughout this 


work^ K will always denote a generic constant ^ich may or may 
not be same* On occasion/ we will show that a constant depend! 
on certain parameters while being independent of others* 
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1^3 Math emat i cal P rel Imln arl es * 

Consider the following parabolic ecjuation in a fixed 
domain I x 

Cl*3*l) -u^. + Lu = f, (x#t) 6 I X (O, t] , 

with either Dirichlet botjndary conditions 

Cl*3»2) u(o/t) = u(l,t) =0, t > 0 

or mixed boundary conditions 

(l*3«2') u (o/t) = u(l,t) =0, t > 0 

and appropriate initial condition, 

where L is a second order (possibly nonlinear) elliptic operai 
Weak formulation . Let ^^(l) be the space of all factions 
<p e H^(l) satisfying either 

<P(o) = Cl) =0 

or 

<P^(o) « <P(l) = O, 

Now the weak solution u of Cl»3<il) and (l*3»2) or (l»3*2'') is 
defined as follows : Find a function u ; -*■ H^( l) such 

that 

u^. e H^( l) in the sense of distribution for each t > 
u satisfies the initial condition. 


and 
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(l*3«3) ^ ® H^(l)/ t. > 0* 

Let us now turn to the formulation of Galerlcin method 
for approximating the solution u of (l*3»3)« 

o 

Galerkin procedure # Let denote a finite dimensional suhsps' 
o 

of Than a continuous time Galerkin approximation of th 

Vi o 

solution u is defined by a function u^ : [O/T] such that 

(1.3*4) + (Lu^ X^) * (f, ^eSj^^0<t<T 

with a suitably chosen initial function u (x/0)« 

Let i = 1/2, •••/N be a basis for the finite 
dimensional space Then the Galerkin approximation u^ e 

can be written as 

( 1 * 3 * 5 ) u^ = L a.(t)vj(x)* 

i=l ^ ^ 

Vi 

Setting the above ej^ression of u and taking X a Vj, j = 1/ 2, 

in ( 1 * 3 * 4 )/ we obtain an initial value problem for the system 

of N ordinary differential equations * These ordinary differen 

tial equations in s are nonlinear if atleast one of L or f 

ll 

is nonlinear in u* Thus, the determination of u amounts to 
solving the system of N differential equations in a^^s with an 
approximate initial condition* This set of equations can be 
solved by the variety of numerical methods available for 
solution of a system of ODE, one of which is to discretize it 
in time and solve a system of difference equations* 


In the 


8 


sequel^ we shall call this the fully discrete approximation « 


ikaproximating finite dimensional spaces* Generally the finite 
... o 

dxraensi-onal space considered to approximate u in the 

Galerkin procedure^ belongs to a regular family and 

satisfies certain inverse properties* We recall her® the 

IC 

definitions of families as described in Babuska and Aziz 

[ 5 J , Oden and Reddy [ 56 ] » 


Definition 1 » 3 * 1 * A class of finite dimensicnal subspaces 
^h* ^ ^ referred to as an family if and only if 

the following conditions are satisfied t 


( i) c r+1 > k > 0* 

( ii) a t where be the space of polynomials of 

degree < r^ for regular refinatients of a finite element 
mesh over the prescribed domain* 

(iii) For each fixed he (0, l] and for each v e there 
exists a constant K^/ independent of h and v such 
that for 0 < j < min (m,k) and r > 0 


( 1 • 3 *6 ) inf 1 I V-X.I 1 j 1 I V i 1 

xes^'^ 

where a = min(r+l~j#m-j ) • 

For exarrple of such spaces when k = 2 # one can refer to 

CiarlGt[l 3 j, Oden and Reddy [56]# Douglas et al« [20] and 

Ic ^ 2 

12 J* space contained in H will be 

called spaces* 
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De finition l»3*2* (Inverse property)# A family of finite 
dimensional subspaces is said to possess the inverse 

property if there exists a positive constant such that for 
all j < Ic, 

(1*3*7) 1 1X1 ixl 1 j 4 for X e # 


In the one dimensional case, we shall sometimes refer 
to the following as the inverse property (it will be clear from 
the cases) 


(l»3*8) 


J 1X1 


W' 


3,00 5 


K_ h 


ri/2, 


lx: 


0 < j < X, 


for every X s 

The assTimption that such a property holds for all 
elements in can be found in Nitsche [4?] and Ciarlet [l3]* 

From time to time, we assume one of the following inverse 

o ^ 

properties for the space which belongs to a regular ' 

family 


(l.3#9) - ^o 

and 

(1.3*10) IIX^JI, for X e s^* 

Iti 

We are interested in finding the error U'-u in terms 

ll 

of h» This estimate measures the closeness of u to u in 
suj-table norms* If we get an error estimate for the Galerkin 
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approximation in norm with the same power of h that is 

present in the approximation property of the finite element 
o 

space we call the error estimate an optimal estimate* 

In the following we recall a well-known Theorem for 
establishing the existence and uniqueness for operator equation 
s at i s f y ing co ere iv ity [l 3 J • 

Theorem 1»3*3* (Laxi-Milgram Theorem)* Let V be a Hilbert space 
with norm i I *| i and A(«, ,) : V x V -* R be a bilinear form with 
the following properties : 

/ 

(i) boundedness of A(*, ».) : lA(u,v) 1 < Ml lull llvll 

(ii) V-coercivity of A( ,, *) : A(u^u) >a Mull^, a>0jU,vev* 

Then for any continuous linear functional f on there exists 
a unique solution u^ in V such that 

A(u^,v) * f(v), for all v e V* 

We now state some results concerning the regularity 
of solution to a general class of elliptic boundary value 
problems of the form, (see Oden and Reddy [ 56]) s 


Lu = f in 1/ 

( l*3*ll) 

B^u =0 on 31, 


j 4 j 

where Lu = S (—1.) D (a.. D-^u), 

i,j=0 

= S b. . D^U , O < < 2nwi 

and f is a oiven function* 

Ik. 
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Let L"' be the adjoint of L and A(u^v) bo the 
bilinear form associated with the elliptic operator L-* Consider 
D(l) = { <P e H^(l) : Bj^<P = 0, X e 31} 

For V s d(l), 

A(u^v) = (Lu,v) s= (Uj.L*v)«i 

Assume further that Ci) A(u^v) is bounded, that is 
lA(u,v)l < Ml lul \^\ Ivi 

(ii) A(«,,) is strongly IiP-coercive that is 

A(u/u) a > 0 and u,v S 

Now we state- the following regularity ‘rheorer% see Oden , and 
Reddy [ 56 ]/ Necas [ 46 ] and Agnvon [ 2 ] ♦ 

Thaorgn i*3a 4» If f S k > 0 and A( *, •) satisfies the above 

tvTO properties, then the solution u of 

(l»3*12) A(u,v) = (f,v), for every vS D(L) 

belongs to and there is a constant independent of 

but depending on m such that 

(1.3.13) 

Some inequalities * Below, we present without proof inequalities 
vrhich will be of frequent use in error anal^7sis in subsequent 
chapters • 

P roposition l-3*5» If ^ s H^(l) and <p vanishes at any point 
x' s [ 0 , 1 ], then 
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( 1»3»14) I I9>1 1 < 1 I I 

and 

(l»3*15) l<»>(x)l < V2I l<Pl I® for St' e [o,l] * 

This is one dimensional Poincdre inequality and some times known 
as Rayleigh- Ritz inequality/ see Hardy et al# [asj. 

Proposition i.*3#6* Let then 

(1<.3.16) 11911 „ <11911 + V2 1 1 91 f^^^l 19^1 

l“( I) “ 

Lf>t us recall here the Sobolev imbedding theorem, the 
proof of which can be found in Adams [ 1 ]/ Oden and Reddy [56]* 

Proposition t*3»7* Let Q ?oe a bounded and smooth domain in 
R^* Let j and m be non-negafcive integers and let p satisfy 
1 < p < os>. Then for mp > n, an imbedding of the form 

T/^+j/P(Q) ^ w^''^(q), p < q < , 

holds that is there is a constant C = C(m/n/p/q/Q) such that 

11911 <011911^,. #p<q<<» and where 

W^'^(q) and 1 1 • 1 I ^ are the Sobolev space and the norm 

W”'P(Q) 

respectively, for detailed description refer to Adams [ 1 ] • 

In one dimensional case, the following proposition is 


easy to establish ■» 
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Proposition i*3*8* Let (p B then 

(t.3*l7) 1 1^1 I «, < Cl |<pl i • 

L ” 

Further* if (p vanishes atleast at one point y/ e [o*l] * then 

( 1 *3*18 ) I 1^1 I ^ < 1 I 1 . 

L ^ 

We shall also use the following inequalities* For 

more detail, see Hardy et al» [28 ]« 

2 2 

I* ^ ^ 3,,b > 0 and s > 0» 

aP 11 

II* # for all a,b > 0, p > l and p q ~ 1 

( Young' s in eqnality ) • 

In integral form, if f, g are both real valued and f e L^, 
q 

g e L , then 

^ J ‘ ^ ‘ + CT * * ^ i i “ ^ P > i • 

I ^ ^ PH 

III* (generalization of Young's inequality)* 

aP b^ 111 

abc < ^ + -^ + -j-i for all a,b,c P ^ ^ 

Similarly in the integral form, for real valued functions f,g,h 
if f e L^, g e and h e l^, then 


/ fgh <i J 

I P ^ 





l Ihl 


4 . ^ ^ — 

P q r 


1/ p > 1* 


IV* ^ fg < 1 Ifl g foir fBlP/ gBh^ and p q * 

p, > t (Hardy's inequality)* 
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For p = g = 2f this Inequality is known as Schwarz* s inequality* 

We state without proof, the following Lemma* For a 
proofs see Halanay [27, pp* 7 ] * 

Lenroa l*3it9 (Gronwall's Leittna) <» Let f,g and h be piecewise 
continuous non-negative functions defined on an interval 
^ 5 9 being non-decreasing* If for each t e [a-,b] 

t 

f(t) 4- h(t) < g(t) +1 f(t' ) dt' , 

~ a 

then 

f(t) + h(t) < g(t) e 3 <p (t-a)- 

The following is a discrete analogue of Gronwall*s Lemna, see 
Lee [ 39 ]• 

Lenma 1*3*10 (Discrete Gromvall^s Lemma) • Let f(t), g(t) and 
h(t) be non-negative functions, defined on Ta = {t e [o,t] : 
t = jAt, j r= o, l,##.,N/ NAt = T} , g(t) being nondecreasing* If 

t-At 

f(t) + h(t) < g(t) + CAt S fC'c), 

T=0 

where C is a positive constant, then 

f(t) + h(t) < g(t) ejqj (Ct)* 

We close this Section by stating the well-known 
Schauder fixed point theorem* For reference, see 
IstrKtescu[ 32] • 

Let X be a Banach space and Dom(f) =s Done. in of f« 
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Theorem (Schaudar fixed point theorem)* Let C be a 

bounded closed convex nonempty subset of X* Let 

f : C O Dom(f ) - C 

be completely continuous* Then there exists atleast one u s C 

o 

such that 

f(u ) = u • 

' o o 

1 *4 A Brief Survey * 

For the last several years, the free boundary problems, 
particularly Stefan problems have received steady attention on 
account of the emergence of new areas of applications providing 
significant impetus to their study* This is guite evident from 
the fact that during the last decade many international 
conferences were devoted exclusively to this area of mathematics* 
Below/ we present briefly a survey of finite element methods 
for numerical solution of free boundary problems* Several 
detailed survoys on numerical solution of free boundary problems, 
especially Stefan problems are available due to Nitsche [43], 

[ 51 ] , Furzeland [ 25 ], Mayer [ 43 ], Elliot and Ockendon [22 
and Tarzia [59]* However, the presentation and direction in 
the present survey is more geared to the type of equations, we 
are dealing with* 

The techniques, used for numerical solution of the 
free boundary problems can conveniently be grouped into 'fjront 
tracking methods', where the position of the free boundary is 
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predicted along with the solution of the governing equations 
and' fixed domain methods', where the reformulated governing 
equations are solved over a fixed domain and the position of 
the free boundary is determined a posteriori# 

Front tracking methods » The front tracking methods require that 
the phase front evolves smoothly in time and ;^aee,.and cme usuall 
has to draw on some a priori knowledge of the solution based on 
physical model in order to judge whether the front is trackable* 
Of course, in a given application such decisions are usually easy 
to make* Let us now look at some specific front tracking methods 

Bonnerot and Jamet [ 7 J,[ 3 ] proposed a finite element 
front tracking method in which they discretized the domain by 
means of isoparametric space-time finite elements* The free 
boundary is approximated by a polygonal line whose vertices 
coincide v/ith triangiolation nodes* In this method the change 
of the free boundary is tracked continuously,' while the elements 
deform continuously* It app'Sars from the numerical experiments 
that the method is of second order accuracy* Recently t 9 ] , 
they have extend-ed to two dimensional problems, taut the order 
of accuracy is one instead of tw* In an atterrpt to improve 
the accuracy of the method, Bopperot and Jamet [ 10 ] proposed a 
similar space-time finite element method for one dimensional 
Stefan problem with higher degree (biquadratic. Instead of 
bilinear) elements* It provides approximations which are 
continuous in space variable, but admit discontinuities with 


17 


respect to time variable at each time step* The accuracy is 
seen to be of third order# as seen from numerical es^eriments# 
although no mathematical proof to that effect is available* In 
addition to accuracy# this method is specially appropriate for 
the computation of the solution# which admits singularities at 
the initial time or on the boundary* Purther# in their work 
in 1931 they [ 11 ] extend and modify this third order accurate 
front tracking method in order to solve one-dimensional multiphast 
Stefan problems with appearing and disappearing phases* The 
numerical results show that the free boundaries start smoothly 
with a vanishing initial speed which has not been mathematically 
proved* 

An analogous method# suggested by Lynch and O'Neill [ 40 _ 
incorporates a continuous mesh -deforming scheme directly in the 
approximating finite element equations to take into account the 
displacement due to the change of the free boundary* The basic 
difference between this and Bonnerot~Jamet method lies in the 
fact that the present method is based on the Galerkln approxi- 
mation in space and uses finite difference in time#, where as 
the other involves space-time Galerkin approximation* By 
allowing continuous mesh deformation# as dictated by the boundary 
condition# the free boundaries lie on element boundaries* This 
circumvents the difficulties inherent in interpolation of 
parameters and dependent variables across the regions where 
those .quantities change abruptly* 
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Mori r44 ] suggested a luitped mass finite element 
method, based on time dependent basis functions for solving on® 
phase Stefan problem in one space dimension* The stability 
and convergence of the method has been studied without any erxror 
estimates and some numerical exan^les are provided to illustrate 
the method* This idea is carried over to two phase problem and 
to problems in higher space dimension, see Flori [ 45]# Like all 
the above methods, the mathematical results concerning the rate 
of convergence for the present one are still not worked out* 

Alexander et al* [ 3 j .recently applied the moving 
finite elomont method with piecewise linear elements on trxangles 
with moving nodes to the Stefan problem in enthalpy form* 

Rubinsky and Cravalho [57] also presented a one dimensional 
finite element method in which the grid is stationary and a 
moving node that coincides with the •oo sit ion of the change of 
phase interface is tracked continuously in time* 

Jaisuk and Rubinsky [33] developed a numerical method 
involving finite elements for the solution of multidimensional 
heat transfer problem in the presence of phase change* The basic 
idea of the method is as follows : The energy equation in the 
media and the energy balance equation on the change of phase 
interface are talcen as independent governing equations and are 
solved by finite element formulation, yielding the teirperature 
distribution in the media as well as the continuous displacement 
of the interface* 



19 


The applicability of the aforementioned spac&^time 
fi.nite element method of Bonnerot and Jamet is limited by the 
requirement that the jump condition must be of explicit type* 

In more than one space dimension^ their strategy of determining 
the movement of the free boundary meets with difficulties* 
Following their idea of discretization of the weak formulation 
with isoparametric space-time finite elements* Hsien Li [ 30 ] 
has proposed quite recently a method starting from a weak 
enthalpy formulation* It avoids explicit treatment of the junp 
condition on the free boundary so that the computation is 
simplified and the method is more flexible* It can deal with 
conditions more general than those which are ancountered in 
the enthalpy method* Let us remark here that although 
rigorous error estimates have not been done yet, an analysis 
for a one dimensional two phase Stefan problem does indicate 
that the discretization schema with biquadratic elements is 
convergent and unconditionally stable* A second order converg- 
ence can be expected under some smoothness assumptions on the 
solution and on the free boundary* 

Let us now turn to the method of lines for the 
finite element solution' of the free boundary problems* This 
is based on reducing the problem in one space dimension to 
a sequence of ordinary differential equations with free 
boundarias-Finn and Varoglu [24] developed a finite el orient 
technique, based on this to solve a one dimensional single* 
phase Stefan problem* 


It employs a functional \^dth variable 



domain of integration* However^ this work does not contain 
anything dealing with the question of the rate of convergence* 
There is a need for further research in this direction* 

In [6o] Wellford and Ayer discussed a space-time 
discontinuous finite element method^ introduced earlier by 
Wellford and Oden [ 61 ] for a multiphase heat conduction problem* 
The method uses a fixed grid of standard space-time finite 
elements in conjunction with certain special space-time finite 
elements including the free boundaries* The finite elranent 
model on special elements is defined by discontinuous inte30- 
polation* The resulting model is a true free boundary formula- 
tion in that the jurtp in the heat flux and the position of the 
free boundary are dependent variables* Only numerical results 
are discussed xvithout the rate of convergence* 

Fixed domain mvethod* Lack of smoothness, required for a front 
trackina method makes it difficult to apply in certain situat- 
ions • The alternative is to essentially ignore the free 
boundairy position by solving a reformulated problem over a 
fixed domain* 

The most oopular recent fixed domain method uses a 
reformulation of the free boundary problem as a variational 
or quasivciriational inequality* This approach is still under 
active developmesnt* In the following, we give some recent 
publications in this directions* 
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Ichikawa and Kikuchi [31 ] ( [37] are concerned with 

the fonnalation and finite el^nent analysis of one or two 

phase Stefan problems by freezing indej?;, which is obtained by 

special transformation of the temperature field* The resulting 

variational inequality has been discretized into a system of 

linear inequalities/ \^hich can be solved by optimization 

t echniqueslike the proj actional successive over*- relaxation 

( SOR) method* The convergence of this method has been checked 

numerically* In an attempt to present error analysis Oden 

and Kikuchi [ 55 j considered a finite element approximation 

of a one-phase Stefan paroblem encountered in the study of 

thawing of a frozen media* For implicit finite difference in 

time and the linear finite elements in space# a second order 

2 

convergence in L -norm is established* 

Elliot [21 ] using finite element methods has expressed 
the elliptic variational inequalities# resulting from the 
implicit time discretization of the weak formulation of the 
two-phase Stefan problem in terms of a quadratic progranming 
probl?m\/ which can be solved using proj actional SOR or 
conjugate gradient methods* The SOR method is shown to be 
globally convergent* There is a good survey in this regard 
in Elliot and Ockendon [ 22 ]• More abstract theorems on the 
convergence of the finite element solution of parabolic free 
boundary problems can be foxind in Jerome [34] and Gastaldi [ 26]* 
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Another irrportant version of the fixed domain method 
is the so called onthalDy fo3rmulation for the Stefan problem in 
conductive heat transfer with change of phase* Here the 
equation is written as 

(l*4.i) v.CkVu) ~ = q , 

where H(u) is a monotone function of u with a junp discontinuity 
at the phase change temperature (usually taken as u = o)* 

The numerical solution of (l*4*l) is reasonably well understood# 
One can for example, smoothen H(u) and apply standard solution 
techniques for the resulting nonlinear diffusion equation* 
Smearing the jump of H over a non-zero range (-8,c) amounts 
to a physical model where the substance undergoes a phase 
transition in this zone rather than at a fixed terrperatu 3 re 
u = 0 • The free boundary is identified with the level set 
u = 0, available from the solution of The first 

num.erical vXJrk in the two phase Stefan problem based on 
enthalpy formulation of which we are aware was the 
constructive existential analysis of Kamenomotskaja [36]- 
It was based on an explicit finite difference scheme. Another 
work of the same nature is the result of Milton Rose* During 
l960's there appeared other papers modelled on Kamenomotskaja 
viz* Samarskii and Moiseenko, Budak, Soboleva and Uspenskii* 

One must mention here the work of Meyer [42] along with that 
of Solomon and Lazaridis in this connection* For a recent 
survey, see [ 43 ] and [59 ]* 
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^ [ 35 ]/ Jerome and Rose have der5.ved the rate of 

convergence for regularizations of the multidimensional two 

phase Stefan problem with homogeneous Neumann boundary condition 

and used the regularized problems in enthalpy form to define 

a backward difference scheme in time and Galerkin approximation 

in space over a c'^~piecewise linear spline space* They have 

establianed an L -convergence of order '/’s for the e— regular!— 

zation and an L -rate of convergence of order (h /6 + At/Ze) 

for the Galerkin estimates* This leads to the natxiral choices 

of s At ^ and a resulting 0(h^’'^^) L^— rate of 

convergence of the numerical schane for the v/eak solution o.f 

the enthalpy equation* An essentially 0(h) rate is demonstrated 

2 

when e = 0 and ^t h in the Galerkin schema under a bovmded- 
ness hypothesis on the Galerkin approximations* The latter 
result is consistent with the computational experience. 

Yet another approach to ('l»4*l) is to use u(h) in 
the enthalpy formulation leading to equations 

RCP 

( 1 -4-2) ^ = AA(<P) + Q , 

where u = A<p =h*^( 9) and Q = q o H The phase transition 
now takes place at u = O-but the so called mushy zone appears 
at least in the numerical solution of (l*4#2), where both 
phases coexist* Clavaldini [ 14 ] discretized the weak form of 
(l*4*2) by a quadrature rule and solved the resulting problem 
using both explicit and iitplicit finite element sch< 2 mes* 
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ConvergancG proofs for the explicit finite elment schemes were 
given by Ciavaldini and for Implicit schemes by Schafer in 
Hoffman [ 29 , III] • Recently Epperson [ 23 ]# based on the 
weak form of (l. •4»2) developed a finite element method for a 
nonlinear Stefan problem. Semidiscrete error estimates have 
been established assuming only Lioschitz regularity of the 
nonlinearity* The estimates obtained here are of lower order 1/2 
as opposed to 2/3 in [35] » 

The last fixed domain method in i-jhich our interest 
lies can be treated at the same time as a front tracking method 
because the free boundary explicitely enters into the differential 
equ.ation# This is based on a co-ordinate transformation to 
reduce the problem to one in a fixed domain* The penalty one 
pays for it is the introducticxi of a highly nonlinear term into 
the transformed problem* This requires considerable regularity 
of the solution* The idea of co-ordinate transformation to fix 
the free boundary was first proposed by Landau [ 38 ] and first 
used numerically for a finite difference approach by Crank, see 
Ockendon and Hodgkin [ 54 ] • Since its initial introduction, 
the idea has been used widely in differing numerical techniques 
to solve the resulting fixed domain problems* For example 
Budak et al*, Wentcel/ ai^id Uspensky used finite difference 
schemes for nonlinear Stefan problems both for single phase 
and multiphase situations* For references, sec Tarzia [59]* 

Based on this, Nitsche introduced in a series of 
significant papers [ 48 ] - [5 3] the study of error analysis for 
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th 3 finite slcment Galerkin approximations to the single phase 
linear paroholic free boundary problems in one space dimension* 

In [ 48 ] / a semi-'discrete Galerkin method for a single phase 
Stefan problem was proposed and analyzed* The optimal error 
estimates were also derived* In [ 50 ] ~ [ 51 ] / the error analysis 
was further extended to cover the case such as those which 
occur in oxygen diffusion problems* Nitsche [Slj also discussed 
some refinements and fuirther generalizations, especially related 
to the Problem III in Maganes [ 41 ] • 

1 *5 Outline of the Dissertation* 

The layout of the thesis is as follows : In the second 

chapter, we deal with a single phase non-linear Stefan problem 

with Neumann boundary condition and the same equations with 

Dirichlet boundary conditions forms the object of study of the 

third chapter* In the second chapter optimal rates of 

1 2 

convergence in H and H -norms are established for continuous 
time Galerkin approximations, where as in the third chapter 
optimal error estimates are deduced for continuous as x^ll as 
discrete time Galerkin approximations* In these as well as 
the two subsequent chapters, the approach has been to 
straighten the free boundary by a suitable co-ordinate transfer— 
mat ion and then to apoly H -Galerkin procedures for approximating 
the solution* Further, the global existence of the Galerkin 
approximations has been established by fixed point arguments* 
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In -the fourth chapter, a singlephase Stefan problem 
with quasilinear parabolic equation in non— divergence fom is 
considered and a priori estimates in L , H as well as H norma 
are derived* 

Chapter five describes both fully and semidiscrete 
H^~Galerkin approximations for a semilinear Stefan probl®n with 
a nonlinear source term* Unlike in the previous chapters^ this 
problem leads to a system of coupled equations* Therefore, a 

more careful treatment is necessary* Optimal rates of conver- 

2 1 . 1 2 

gence in L , H , W ' and H -norms are established for this 
case as well* 

In the sixth and concluding chapter, an overview of the 
results is given- The problem of sup erconvar gence and various 
modifications of the Galerkin methods for other related problems 
are discussed* In particular, the Galerkin approximations for 
one dimensional ablation problem and porous medium equations 
are investigated* 



CHAPTER 2 


NONLINEAR STEFAN PROBLEM WITH NEUMANN BOUNDARY CONDITION 
2*1 Int rodxict Ion » 

In this chapter we confine our discussion to the 
singlephase nonlinear Stefan problem with Neumann boundary 
condition , Earlier, Nit sche [48 ] considered a linear problem 
and developed the error analysis for the same. The present 
work is an extension in the same direction to a class of 
nonlinear problems. Here we fix the free boundary by a suitable 
transformation, which results in decoupling the equation into 
two main systems, leading to considerable simplification in 
analysis inspite of the presence of the nonlinearity. 

In Section 2, we state the basic problem along with 

regularity conditions. We discuss a co-ordinate transformation 

to fix the domain in Section 3» Section 4 deals with the weak 

formulation and an H^-Galerkin procedure for the transformed 

problem. In Section 5, the Galerkin approximation is analysed 

for the local existence and uniqueness* An auxiliary projection 

is discussed in Section 6 alongwith the related error estimates. 

The main results are presented in Section 7. where a priori 

error estimates for the continuous-time Galerkin approximations 

are established together with the global existence of the appros^ 

imate solution of the original problem. It is shown that the 

1 2 

order of convergence are optimal in H and H —norms . 
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2*2 Statement of the Problgn» 

We consider the following single phase nonlinear 
Stefan problem in one space dimension^ arising in melting of 
solid or freezing of liquid [58] • 

Problem P : Find a pair C(U,s) : U = S = S(t)} 

defined on the closure of the parabolic domain 

Q X (0,T^]^ where Qr={y : 0<y< s(t )} and 0 < T < 
such that it satisfies 

(2.2.1) ~ (a(U)U^)y = 0, (y^-r) e q x (0,Tj 
and initial/ boundary conditions 

(2.2.2) U(y/ 0) = g(y)/ y e I; 

U^^(0,-O = 0 , 

( 2.2.3) " T > 0 

U(S(t),t) = 0 / 

together with the free boxindary relations 

(2.2.4) S^ + a(U)u 1 = 0/ T >0 

^ y’y=S(-c) 

and 

S(0) = 1. 

Here/xve make several assumptions which will be collectl~ 

«N> 

vely called condition 
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(2*2«5) Condition 

( i) The pair {U, S} is a unique solution to (2»2*l) - (2»2#4- 

(ii) For e R, 0 < a < t where a is a positive 

constant* 

(iii) The function a( »), i^^iich depends on U only belongs 
to C (r) with continuous and bounded derivatives* 
I'toreover, there is a constant > 0 such that for 

s R 

3q^ 

(iv) The initial function g is assumed to be sufficiently 
smooth and satisfies the compatibility condition, 

gy(0) = g(i) = O. 

Further, the solution {U, S} satisfies the follov/ing regularity 
assumptions s 

R^* U e L°°(0 ,T^,-H^'^^(q(t))) n l“( 0 ,T^/H^(q(t} )) 

nw^-'^'(o,T 7 H^’'’^(q(t) ) ) n w^'‘”(o,T /h^CqCt))) 

0 o o 

and 

S S W^'“(0,T^) , 

CO 

Let K 2 denote a bound for the norms of the functions in all 

ro 

of the spaces in R 2 » 

2 » 3 Straightening the Free Boundary # 

In this section, we fix the free boundary’’ using Landau- 
type transformation [ 33 ] 
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(2*3*1) y * X S('r)> t > 0# 

Further/ in order to decouple the resulting system we introduce 
an additional transformation in time scale/ given hy 

(2.3*2) 'r(t) =/ S‘^(t') dt'. 

o 

If we write U(y/T) = u(x/t) and S(t) = s(t) in the transformed 
variable, then 

(2-3.3) U^(y,'^) = u^(x,t) IP = u^(x/t) s‘*^(t) 

( 2-3.4) U^yCy*"^) = -^xx 

and 

( 2 . 3 »s) u (v t) — 9^ 4 . .j, 9^ 

^2.3.5; - 3 ^ ^ ^ ^ ar 

* vi^s"'^(t) - u^ys"'^(t) ^ s'^^(t). 

But/ 

(2.3.6) = if ^ 

and from (2.3.3) 

(2-3.7) Uy(S(T)/t) = u^(i/t) s“^(t). 

Using ( 2 . 3 . 3 ) ~ (2-3-7)/ the Problem P gets transformed into 
the Problem P. 

P roblem P ; Find a pair CU/S} satisfying 

(2-3<«8) ~ (a(u)u^)^ = -a(u(l))u^(l)x (x/t) e I x (O/T], 
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with the initial and boundary conditions, 

( 2*3*9) u(x;,0) = g(x)/ x e I; 

(2.3*10) u^(0,t) = u(l,t) = 0, t > 0 

and 

(2*3*ll) ^ a(u(l)) u^(l)s, t > 0 

with s(o) = 1* 

Cl early, t = T corresponds to T T^* Mote that all the 

TO 

regularity assumptions for £U,S} easily caxiry over to 
{u,s} and we call these with the bound -^s u(l) = 0, 

a(u(l)) = a(o), a known constant independent of u* Further/ 
the integral (2*3«2) can be rewritten as 

(2-3*12) ^ = s^(t)/ t > 0 

with T(o) = 0* 

Remark - The transformed Problem P involves a nonlinear 
parabolic problem ( 2*3-3 ) ~ (2-3*lO) in a fixed domain and two 
initial value problems for ordinairy differential aquations 
(2*3*11) and ( 2*3*12) • Thus the computation of u can be 
carried out independent of 's'* 

2*4 Weak Formulation and Galerkin Procedure * 

Consider the space s 

H^( I) = { V e H^( l) : v^(0) v(l) = 0} * 
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Multiplying both the sides of (2*3*3) by -*v^ and integrating 
by parts the first teirm with respect to xy we get 


(2*4*1) ^ ’^x^^ ^^^x'''^xx^' 

V e H^(l) and 0 < t < T* 

If u e l'^( 0,T;S^( I) ) with e b^(O^T; H^( l)) satisfies (2*4*1) 
and initial condition u(xjO) = g(x)# then it is called a weak 
solution of (2*3*3) - (2*3 *10)* 

1 

H -Galerkin procedure * For each of a family of values of h 

/ T ® 02/ \ 

in (0,1J, let ^ H*(I) be a finite dimensional subspaca 
satisfying both approximation and inverse properties ( 1*3*6)/ 
k =5 2 and ( 1*3*9) respectively* For exanple of such spaces 
one can refer to Nitsche [ 47] • 

Now we call u^ : (O/T] -» an H-Galerkin approximation 
of u, if it satisfies 


(2*4*2) (uJ^,X^)+((a(u^)iaJ)^,y^^)= a(0)uJj(l)(xuJ/X^)/ Xe 


and the initial condition 


(2*4*3) u^(x/0) = Q^g(x)/ 

o 

where is an appropriate projection of u onto at t = 0/ 
to be defined later* 

Further, the Galerkin approximations Sj^ of s in (2*3*ll) and 
of T in (2*3*12) are given by 
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( 2*4 *4) 
and 


( 2*4,5 ) 


* 


~a(0) u^(l)sj^ 


f 


with SS 'X 


dt 


s^(t ) 


with = 0, 


2»5 Existence and Uniqueness of the Galerkin Anproximation » 

Now we show the existence and uniquariess of the 

ll 

Galerhin approximation u satisfying (2'»4»2) — (2»4*3)* Since 
o 

is finite dimensional, the equation (2'‘4*2) results in a 
system of ordinar^r differential equations in the temooral variable 

t in explicit form with differentiable right hand side* There— 

Itl 

fore, for a given u (x,o) there exists a unique solution to 
(2*4*2) in a certain interval (0,tj^)* Ws shall now show that 
tj^ does not depend on but may depend on the data g» 

Following Nitsche[48 ] , we have 

Theorem 2*5 *1, Assume that the projection is bounded in 

H^-norm* Then there is a t > 0 independent of the choice of 

o h 

the approximation space such that u is uniquely defined 
for t e (O,?]* 

Proof * In order to establish the existence of a t independent 
of S^, choose t e (0,tj^)* Setting X = u in the equation 
( 2*4*2), we get 




But 
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and 


(2*5.3) a(0)u^(l)(xu^,u ) < ‘ eoHu^} I Hu^^H 


Since ■>J-^(0) = O, then we get by ( 1*3.15) 


( 2.5.4) I lujjl \ ^ <'fl \ \v^\ 1 * 


By ( 2 . 5 . 1 ) - ( 2 . 5 . 4 )/ we now have 


dt 


I lujjl I ^ t 2al lu^l 1 ^ < 4N/’2K^ 1 Iv^l 1^/^ 1 1 


Applying Young's inegiaality (Inequality II of Chapter l) with 

h 3/ 2 

a suitable choice of e and q/ where b = el lu I I so that 

XX 

V» o 

2a 1 1 = we obtain 

in • ii:^l I ^ < K2 (Kj^) llujjll^. 

h 2 

Therefore^ X = 1 lu^l 1 obeys the differential inequality 

it - ^3^^^ t > O 


with 


Mo) = I {Qh^x* ^ ^^ 4 ' ‘^x’ ' 


Integrating with respect to t fr6m (0/t), we get 

X^(0) 


- ^ x^(t) < 


l-.2K2t \^(0) 


Prom which it is clear that I lu 1 1 is uniformly bounded in any 

iX 


compact interval of [ 0/ 




■^) which is independent of S^i 
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Here the constant « K^Ck^/K^)* Hence is defined in the 
above interval • 

We have seen earlier that for a given choice of a base^ 
the aoproximation u is chosen uniquely* However, the following 
uniqueness theorem shows that the approximate solution is 
ind^endent of the choice of bases* 

Theorem 2*5*2 (Un3.queness of the approximate solution) • For 
any K > O, there is atraost one solution u^ S of (2*4»2) — 

( 2*4*3) lying inside a ball where 

Bk = {w( -,t) e H^Cl) : 1 Iwl I ~ * 

Proof * Let u^ and u^ both belonging to be two solutions of 
( 2*4*2)# satisfying the same initial condition* Then the 
difference <jp = '^.^“'’^2 S9''tisfies. 


( 2*5*5) 






^xx> 




( 2*5*6 ) qi(x,0) a 0- 


h h ^3 a(u^-§P) 

Since a(u^)~a(u 2 ) = Vd| , the equation (2*5*5) 


can be rewritten as 


(2.5.7) ('P^.X^)+(|^(aCu^)*^),X^)=(|j(/ ^ ^ '»<«<x>'Xxx> 
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+ a(0) tu^ 


Since 9 B we can set X =i<P ^ obtaining 


(2.5«) -j 35 IKP^II +(5;(a(ui)V'‘'xx’“'|5<' 

o 


^ aa(u^~g(p) 






^1 "^2* 


But as In (2*5*2) 


(2.5.9) (|^ >.an<f_l|tK,llu, „I1 „ll»„ll ll<P._.ll. 


XX" l'’“l^x' 


Using I 1X^1 I < I iX^^i I and I Ixl I ^ < I 1X^1 I for X e S , it i; 

L L 

easily seen that 


( 2*5*10) 


151 < l<Vx“ix'‘*’xx>> + ''^<l’'^ixx'''xx>' 


+ 1 (/ 


1 3\Cu5;-|(P> ^ 


o 9 u' 


^"l,x-S»x>''“2,x “^xx^' 




+ 1 l<P,,l I I KPI I J lu, ^1 1 1 la II 


where 


< lCg(K^,K) IK^^l 

1 aa(u^-g<p) 

= / rrr <35 . 




Further, we have 
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(2.3.11) II I < la(0)u!i' (l)(x(p )l+la(0)<p (iXxu!} 

^ J-/X x: XX X XX 

"< x '' 

Lf 

< Kg(K^,K){||<P^ll + |<P^(l,t)l} IKP^II. 

From (2»S*8) (2*=;»ll)* we obtain 

(2.5.12) I if 3Kg(K^,K){)|(|)^lI 1 l<P^I I 

< 3Kg(K3^,K){||<p^ll 
+ /2 * 

Using the Inequality I of Chapter 1 for the first term on the 

right hand side of (2 *5 *12) and Young* s inequality (Inequality 

II, Chapter l) for the second term with appropriate choices of 

2 

e, p and q, we get the terms with 11 (p 11 cancelled from both 

zxx 

the sides and get 

(2.5.13) 1 Icp^l 1^ < K^(Kj^,k) lltp^ll^. 

As <p(x,o) s= 0 and consequently <P„(x/0) = 0, we have from 

X 

the inequality (2.5.13) 

(p^i X/t ) S 0. 

The later shows that <P(x,t) is independent of x and vanishes 
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at X = 0, which shows that <p(x/t) s 0 a»e* This completes 
the proof of the uniquaness of the approximate solution- 

2*6 Some Approximation Lemmas * 

Let us now introduce the following forms : 


(2*6-i) B(u/V/w,X) a(u(l/t) ){ v^(lyt)(xw^)^^) 

+ w^(l#t) (xv^, 

and 

(2*6 *2) A(u,v?w) = (a(u)v^),w^ ) - 2 B(u,u|v,w)* 

The following Leitma shows the 3x>undedness of A and establishes 
0 o 2/ \ 

a Garding~type inequality for A in H (!>* 

Lemma 2«6»l(i) (Boundedness)* A(u;v,w) is bounded in S^Cl) 
t hat is 

(2*6*3) lA(urv,w)l < Ml lv„„l 1 llw„„ll* for u e Vi^*“*(l) and 
o 2 

V/W e H (l), whore M may depend on I lu I I — • 

^ L 

(ii) (Garding-type inequality)* There exists a constant 
a > 0 and P such that 


( 2*6 *4 ) A(u?V/v) > al I v I 1 ^ - P 1 1 v I I for u 6 ( l) and 

o o 

V e H^(l)* Here p may depend on 1 lu 1 I 

^ L 

P roo f » From (2*6*2)/ we have 

(2*6*5) lA(u/V/w)l< I (~( a(u)v^),w^) 1 + 21 B(u/Ujv,w) r 
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< 1 ( a( u ) ) 1 + 1 ( ) I + 2 1 B( u, uf V, w) I • 

But^ using lv(l)| < llv^^ll, for v e H^(l) 

(2*6*6) 21 bCu^u/v, w) 1 < la(u(l) ) l{ lu (l )(xv„/W^ )l 

+ IVx(lUxn^,w^)l) 

<K^{|U^(1)1 iiv^M iiw^ii+iv^(i)i iiu^ii myi) 

Ij 

L 

Substituting (2*6*6) in (2*6*5)^ using (2*2*5) condition ^(ii)/ 
(iii) as well as the Poincare's inequality (l*3*14) llv 1 l<I Iv llj 

X XX. 

O n 

for V e H (l), we get the estimate (2*6*3)* 

For the nroof of (ii)/ consider^ 

(2*6*7) A(u;v,v) 

> (^(^^'^xx''^xx^''^ l(a^(u)u^v^^v^^) l + 2lB(u,u#v,v) 1} 

> allv^ll^-K^CllUxll «,! Iv^I l-fllu^l 1 Jlv^(l)l 

L L 

+ I Iv^l 1 )} I I 1 

1 

> al I 1 lu^l 1 J 21 I v^l Itvr2! I v^l Ilv^^l 1 

L 

Apply Young's inequality for the second term in the bracket 
of right hand side of (2*6*7)/ we get 

(2*6*8) >(2Kj^llu^ll 

L 

< |(0llv^^ll^/^)'^/^ + |(K^ ^>^2llv^l |V2)^ 

" ^ L 
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L 

Further, applying the trivial Inequality I of Chsipter 1 for the 
first term in the bracket of (2*6*7), we get 

4K^‘ 

(2-6.9) 2K^llu^ll ^llv^ll llv^^lKellv^^llV- 

L L 

Now choosing ^ - e and e = a/4 and applying (2*6»s) ~ (2-6.9) 

in { 2*6*7) f we obtain 


A(u,v,v)> |l Iv^l I ^-K^> lu^l I lu^l 1 " I I v^n 

L L (X 


2 (16 . „2, , , , 2 27 

L (X 

r2, ,,2 ~~l(3_6-f27I<'^ 

then the required estimate (2-6-4) follows 


2 

I 

2 -*'2- 


If we identify a/2 = a and K£1 lu^l I a"'-^( l6-f27K£l lu^l a‘'"^) = P, 


Lf=t Ap(u;v,\'^') = A(u7v,w) + Tben Ap is 

o o 

strongly coercive in H^(l),that is 


(2-6-10) A.-(uyv, v) > a I Iv H for u s and v e H^'- 

^ ~ . XX 

For t e [o,T], we define u(x,t) e as the Galerkin aporoxima~ 
tion with respect to the form Ap : 

(2-6-11) Ap(uju~u,X) =0, X 6 Sj^» 

The existence and uniqueness of u can be easily proved following 

the Lax-'Milgram Theorem (Theorem 1-3-3), since the bilinear form 

o 

Ap satisfies boundedness and coercive property in and 
F(X) s= Ap(u/u,x) is linear in X and continuous- Here the 

ISJ 

continuity of F(X) follows from the condition R^- 
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Below/ we derive the estimates for V s u-*u and its 
temporal derivative 

L emma 2«6*2# Let T> = u-u be defined by (2*6*ll)# Then, for 
t e [o,t], th=^re exists a constant K- » K (a/P,M,K ,K^;K-) 

o o OX/ 

such that 


( 2*6 * 12 ) I 1 < Jul , 2 < m < 

Proof* Take X - V (u-v), for v e in (2*6»ll) < By the 

O 2 N 

coercive property (2*6 *10) of in H (I), we get 
^ ^ ^ ~ Ap(uf'n, T?) ss A (u/ri,u~v) 


< M iin^n 

Applying approximation property 
required estimate (2*6*ii)» 


inf llu-vll 
veSj^ 

(1.3*6) for S 


2 * 


h' 


we obtain the 


We shall also require the estimates for the projection 
2 1 

error in L and H -norms* In the following, we use in 
a generic sense* 

Lemma 2*6*3* There exists a constant as in Lemma 2*6*2 
such that 

(2*6*13) 1 im J < Kg h^ 1 1 7^1 1 • 

In particular, 

( 2#6*14) 1 ini 1 < Kg h™ 1 lul , 2 < m < r+1* 
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Proof * The proof will follow from the standard Nitsche-Aubin^s 
duality arguments with suitable modifications# Consider one 
dimensional formal elliptic operator given by 


( 2 # 6 * 15 ) L“(u)w = 1^ (a(u)w^) t a(u(l|t))u^( l,t)( - Pw* 

Let ^ e L^(l) and define (p e h'^(i) H h^(i) by 


( 2*6 # 16 ) 


l'^’Cu) (P^ =1! , X e I ; 


( 2*6.17) 


=! (xu /<P )? 

xxl . x' XX ' 

X=1 


( 2 * 6 . 18 ) 


<P. 


XXX 


= 0 # 


x=0 


Then, for Ve H^(l) 


(2*6*19) (VfL'' <P^^) = (n, •~(a(u)<))y^^))'’r(r),a(u(l))u^(l/t)(x;f'^^) 


XXX 


X 


XX 






~ a(u(l))u^(l)(x ^/‘Pxx^ 


+ a(u(l))u^(l)x)1 <P^]^ + P(V*x> - P'W'xlSo 


X=1 


(Using integration by parts)* 

The second, fourth and sixth term vanish because of (2*6*18) 

, O 9 » 

and T) as well as <P s H (l)» Integrating by parts, the first 
term on the right hand side of (2*6.19) with respect to x, 
and using (2*6*17)/ we get 
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(2.6.20) (P^) = ((aCu)Ji^>^.(p^> - a(u) 


- a(u(l))u^(l) (x 

= (^a(u) - a(u(l)) Tj^(l) <P^3^Cl) 

- a(u(l))u^(l) (xT)^, (p^^) + 

= Ap(uyT), <p)* 

Thus defining D(i;'*' ) as follows ; 

(2*6.21) D(l'-’) = {<P e H^(I) n H^(I) 5«’xx^l/t)=(xa^/<Pxx^ 


(D (0,t) 
^XXX ' 


0}, 


we have a solution cp e D(L'' for each ^6 L (l) satisfying 
( 2*6.16). This follows from the positivity of the operator Ap. 
Here L* is a bounded map on DCl**" ) onto L *( l) with 

C”^l 1 < I I (L^)“^^l i 2 < !«• I • 

Thus^ we have 


( 2 . 6 - 22 ) 


1 1 g)n^ < li-l I, 


where C depends on I lul 1 - • 

o vr' 


Multiplying both the sides of (2.6*16) by Vf we get 


(T)^i) = Ap(uj“r7,^) 

= Ap(u/‘»7,(p-X)., for X s 

<M I I ^ l)<P~Xll2* 

Xegj^ 
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Here we have used (2*6 *20) j (2«6-»ll) and (2»6«3)* 

o 

By approximation property ( 1*3*6), k = 2 for and the regula- 
rity property (2*6#22)> we obtain# 

(2.6*23) (T1, '!) < Mj 1 h^McpIl^ 

< MKoCoh^l 1 1 IS U . 

Our estimate (2*6 .13) for liTjl) follows directly from C 2.6 h* 23) * 
To obtain (2*6.14), use the estimate (2*6.12) in (2.6.13). 

Lemma 2.6*4* There is a constant as in Lemma (2.6.2) such 
that 

(2*6.24) n-ni < Kg h! U 

In particular, 

(2*6.25) I 1-171 < Kq h^^ 1 lul Ijjj, 21m<rfl. 

Proof. The result ( 2*6.24) follows as a conseciuence of inter- 
polation inequality 

I 1771 1^ < 1 Ihl 1^/^ 1 lT?l using I l^Tl 1 < 1 177^1 1 and 

1 I I < II n^x* * '* ^ ® " 

Applying the estimate (2*6*12) (2*6*24), we 

get the rec[uired estimate for (2*6 *25). 

In the following Lemma, we obtain an estimate for 7)^ at x = 1 < 
‘Lomma 2*6 *5'. There exists a constant (o^,,KQ*rKj^.#I4,Kg ) 


■such that 
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( 2,6.26) Ii 7 ^(l,t)l < Kg I )uj , 2 < m < r+1. 

P roof . Define an auxiliary function <P e l) D H^‘(l) 
satisfying 

L"'(u) = 0, X e I# 





XX 


= 1# 


(p. 


XXX 


= 0, 


X=0 


Multiplying both the sides of the first equation by n and 
integrating by parts, we get 


la(u(l)) T?^(l,t) I = IAp(u/i7,<p) I 

s |Ap(u;T?,<p-x) I , X e 

inf ll<p-.X!l 2 

Xeih 

Using the estimate (2*6.12) and a(u(l)) > a, we get the desired 
result . 

We now wish to find out the bounds for in different norms. 

Lemma 2.6*6. For t s [o,t], there exists a constant 
= K^Q(a,P,M,K^,K^,K 2 ,Kg/Kg) such that 


( 2.6.27) i 2<m< r+1. 

Proof . Differentiating (2.6.11) with re:^ect to ^t''. 


we oet 
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(2.6*^) Ap(u/T)^,v) » ^(a(u)))^v^^)+2B(u/ii^/T?,v)* 

o 

Let u'' e satisfy 

(2*6*29) Ap(u;u^''U''‘yv) =0, V e 

Choosing v « (xi|.-u )-.(u^-.x) in (2*6*29), for Xe and applying 
boundedness and coercive property for we get 

a I i ^ ~ A p( uf u^~u*, u^.-u* ) 

= Ap(u;u^-u*, u^-X) 

< Mj 1 (u^-u*)^^! I inf llu^-XIlj* 

xeSh 

Thoraforo, by the approximation property (1*3*6), X = 2 for the 
o 

space 

( 2.6*30) I I I < a~^ MK^l lu^l h^^, 2 < m < rfl* 

Subtracting (2*6*29) from (2*6*ffi) and choosing v = u'**’~u^, we 
obtain 


a I iu"''~u>^. I i < Ap(u 7 u''''''U^/U'' -u^) 

Xj 

Hence, from (2*6*12), (2*6>24) and (2*6*26) 

(2.6.31) llu'f-Stll < S"\Kj(Kgnull^h'^^+K^llull^h'*-l 

+ Kgllul 1„ 
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From ( 2 * 6 * 30 ) and (2*6*31), the result follows* 

2 

The following Lerttna gives the estimates of in L —norm* 

We use in the following K^q as a generic constant* 

Lemma 2*6*7* There is a constant as in Lemma 2*6*6 such 
that 

( 2*6*32) 1 \\\ 1 < h^( I luj.) I^+l lul 1^)^ for 4 < m < r+1* 

Proof * The proof makes use of duality argument in a modified 
form. 

Let 6 D(l‘*') be the solution of 

^ / for W e l^(I)* 

Clearly, by regularity property (2*6*22) H^PU^ < C^I1H;II* 
Multiplying both the sides by we get as in ( 2*6 *28 ) 

( 2 - 6 . 33 ) 

= Ap(u,T^.,<p-)^-(|^ (aCu)^ Tl^, 

+ 2B(u,Uj./TJ,X) 

-Ap (u;ti^,(p-x) + (§ 3 ^ (a(u)^ 

- 2 B(u,Ut/^/<P"X)-(| 3 j (a(u)^ T)^), <P^)+2B{n,n^fV,<P 

= ^2 ^3 ^4 ^ 5 ' 

The I 2 and terms are bounded by 

( 2 . 6 . 34 ) Kj^KjCim^ll + inf I l<p-Xl 1^. 

Xeij, 
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(2*6.35) < ^^2* *\xx" nfl>-Xn 2 * 

Xeij^ 

5‘or integrating by parts we get 

I 4 = 

and (p e D(L*) ► 

Therefore, 

(2*.6*36) II 4 I < K^V'^'^xxx^' + “‘Pxxxx"^ ^ ^ 

< K^K^ I 1 <p| I 4 I IT?| 1 . 

Analogously, we obtain 

( 2 . 6 . 37 ) II 5 I <K^K 2 (IIt? 1 I + lr?^(l,t)l) 1 1 (Pi 1 4 . 

Combining ( 2*6.33) - ( 2*6.37) and applying the regularity 
property ( 2 . 6 * 22 ) and the estimates in (2*6.14) and (2*6*25) 

( 2*6 .27) /we get 

( V»> i 1 ^, ^ 

Kjoh'^^ I l„-)-l lul I J I ICfl l^+KjKjdCg I lul l^h"VKgl lul l„h2("*-2)) 

ll»ll4 

S I lul 1^+1 lut 1 1 „)h"'+Kg I lul ) n$| I 

Therefore, for 2m~4 > m that is r+1 > m > 4 
(2*6.38) (n^,^) < K^q h'^dlull^ + llUtllrn^ IISII. 



49 


and th«9 result follows directly from (2*6*3)« 

R Q.nia.rk » The estimate for 1 I can be proved for m > 3 that 

is r > 2, if we consider ^ = ^txx ^ proof of the Lemma 2*6.7 
or by interpolation similar to the Lemma 2»6*4 between llTt 11 
and 1 |T7^| I , 

2 ■* 7 A Priori Error Estimates for Continuous Time Galerkin 
Approxjunation * ” 

In this Section^ we consider the global existence of 
u and a priori estimates for u~u^. Denote e sa u—u^^ h = u— u 
and I — u^-U/ then e = T)— Let u^(x/0) =s Qj^g be defincad as 
the Galerkin solution of u(x/0) with respect to the form A p : 

Q 

(2»7«l) Ap(u;u~u ,X) ssOattssO, XeSj^ i*e* A^C g-Qj.jg^ X) = 0 » 

Since g is sufficiently smooth/ the projection Qj^g is defined 
uniquely by Lax-Milgram Theorem (Theorem 1.3*3) » Clearly, 
u^(x/0) = u(x/0) • 

Below, we rewrite the Galerkin approximation equation 
in a suitable form. For this purpose, we have frcxn equation 
(2.4*1) and (2.6.11) 

(2.7.2) (u^.^/X^)tAp(u/U/ X) =-'(t 1 j.^/X 3 ^)-B(u,u|U/X)-!-P(u^,X^), Xs 
Subtracting (2*7*2) from (2*4.2) and setting t ~ u ~U/ we get 
( 2.7*3) ( )-Ap(u;U/X)5(T)_j.^X^)-*B(u^,u^/u^/X) 

+ B(u/U7U/X)+ ^h* 
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But, 

(2*7.4) Ap(u^/u^, X) - Ap(u/u,X ) =Ap(u/S,X)-i-C|^([a(u^)-a(u)] u^), 

X„„ )+ 2 B'( u, u/u^, X ) - 2 B ( u^ /u^/ X ) 0 


Further, 


( 2.7.5) 1 ^ ( [a(u^)~a(u)] u^) = ** I3J 


= -; -*§e ) d| e u^ 

o 8 u^ ^ ^ 


«« h h 

+ a eu+a eu > 

^ U X X U XX * 


where 


CM 




3 a ( u— ^ e ) 
3 u 


d|. 


Substituting the values (2*7*4) *- (2.7*5) in {l»l *z) t obtain 

1*1 

(2.7*6) ( Stx^X^)+Ap(u}i:,X) « (n^^,X^)+B(u,ufu,X)-B(u ,u /u , x) 

- P'V’^x’ + 

+ (; ^ ® ''x'’^xx' 

o 

- 2B(u,Ufu^,X) + 2B(u\u^/u^,X) * 

Now, for u and. u^ s we have a(u(l)) ssa(u (l)) s3a(0) and 

(2.7*7) B(u,u/u, X) - 2B(u,u;u^,X) + B ( u^, u^/u^, X ) 

ll 

ss B(u,u^e,K) - B(u,u 7e,x) 

= B(u,e7e, X). 
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Prom (2*7*6) and (2*7*7)/ we get 

( 2*7.8 ) ( ^ "• V^x^ 

+ a(0)e^(l/t)(x Tj^/Xj^^) ~ a(0)e^(l,t) 

+ (/ )^) 

■KJ 

o 

+ a d5 V>7-2)(ax+"x). 

o 3 U ^ 

X 2 

- (/ I d| 

o 3u 

This equation can be interpreted as a.n implicit equation for 
given in terms of t Integrating by parts the first 

and third terms in the right hand side of (2^>7<»8) and r^lscJng 
e by E, for some function E(x/t) e H^(l) for each t/ we obtain 

(2*7*9) ( ^.3j,/\)+Ap(ujS/X) = 

+ a(0)E^(l/t) (x T)^/X^^)-a(0)E^(l,t)(xi;^/>^^) 

+ () d5 (''x"*x^''^x+V' 

o 

+ (/ as (n-s) (E^x+V' Xxx> 

+ (/ 5!|('^-5e) u^(T,-i)(E^+u^), x^) 

o 9u 

- (; 6 « (-ix-Sx^^^V^x’' >^x>* 

o 3 u 
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This is a linear equation in i for each E with 

( 2*7 *10) ^(x,0) =0, 

Hence, there exists a unique solution $ of ( 2*7«9)—( 2*7*10) 

Q 

for any given E e H (l)r Thus, we have defined hy this process 
an operator t on H^( l) through the equation (2»7»9) and 

(2*7*10)* Clearly ^E = TT'-e* 

Vi. 

Now/ thR problem of existence of u reduces to finding a. fixed 
point for the operator e s= 11 ^ otherwords/ since e ?= 

and t depends on e also depends on Em We show tnat there 
is an E such that e = 

Lemma 2*7*1 • For a given E(x,t) s we have the following 

estimate for t in (2*7*9) 


( 2 - 7 . 11 ) Il^jloo 2 



(L^(I)) 



< K. .,£ 1 177. 1 1 ,5 0 +( 1 

~ 11 t 

1 |E| 1 

1 „ . (I-HIEII 

l”(h^) l*"(h2) 

1 I77„1 1 1 2 +1 1 7)1 I } 

^ L^(ir) 

ejq) 



(I'i-l lEl I „ 2 ) + P}], 

where K. . are constants depending on and K^* 

Ju 

Proof . Setting X = i in (2*7*9) and using (2-2*5) condition 
and R 2 we obtain 
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^ Ie ll2^1l\alU^^ll^<lfT)^l| lU^II+pIl tj^ll2-5.pi|T>ii 11^^11 
+ Kj^lE^Cl)! 117^11 I U^l )+K3^1E^(3.)1 llg^ll IU3^35.n 

+ K^( 1 |T>^I 1+1 1 1 )1 lE^l I ^11 ^^^1 \+\K^{ I 1T)^1 1 + 1 1 I ) N 

+K^( I ini l^„+l 1 1 \ i^ooK I lE^^l l+K^)! I ^^1 I+K^^K^C llnl 1+1 I U 1) 

c I lE^i l^„+Kj)l 1 I+K^c 1 ITl^l i+n t^i 1)1 lEi I^J I lE^i |^„+K^)|| 5^11. 

Noting Me II ^ < I I El I - , for E e H^(l) and using 
L “ 

12 2 

ab ^ a '+ s b , a,b > 0, e > 0 to the right hand esq^ression 
with b = I I 1^ we get 

(2-7-12)^ I I ^1 I ^+2al I 1^ 

< 406 11^^1 |^+ I llTi^.11^ + (1 IEII^ 

+ 1 lEI 1 ^(l+l lEl l^)+l) 1 in^M ^+K^3(K^,K^,e )1 inl 1^ 

+ (II El 1^(1+ 1 I El 1^) + P) il 

Here we have used Propositions 1*3*5 and 1*3*8 to estimate ll<pll 

and 1 I cpl 1 oo by I |(,o 1 1 , since <p e I) • 

L ^ 

Choose e so that the term 1 I ^1 I cancols from both sides of the 
equation (2*7*12)* Now, applying Gronwall's Lemraa (Lemma 
1*3*9) and talcing sup remum over the interval [o,JD], we get 
the required estimate (2*7*ll)* 
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The following theorem shows the existence of a unique 
solution satisfying e(E) = E* Purthor# an order of magnitude 
is obtained for the same* 

Thc?orem 2*7«2# Suppose u satisfies the regularity condition 
R 2 * 'I'hen, for sufficiently small h there exists exactly one 
solution u in the neighbourhood of u such that 


( 2*7*13) I lei 1 ^ „ < K. . 4 < m < rfl, 

L (H^) - 14 

where, is a constant depending on K^, 

P roof * in the estimate (2*7*11), using (2*6*14)/ 
( 2 * 6 * 32 ), we get for m > 4 


and 



( 2 * 6 * 25 ) and 


( 2.7.14) , I ^1 +1 ,u, 1^2, 

Now, using the inverse property (1*3*9), for the sp3-ce and 
I U 3 J.I I < 1 I ^xx* ^ ® ^h 
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l2(h"') ^ l2(h^) l2(h«') l”(h2) 


(l+l I El I 


l)!]e=P U+llEII^„^^j^Up,] 


In the bail 1 1 El I ^ < 1, we have 

L (H^) ~ 


( 2 * 7 . 15 ) Dell „ 2 ^ ^ n lul I „ ^ - „ 

l'»(h 2 ) “ 15 l 2 (h^) 


+{ I !ul I 


+ {K3^2T(2*p)}]. 

For h sufficiently small^ 

Hell < 1. 

L (H^) 

Thus, the map E -* e is mapping the unit sphere into itself, for 
sufficiently small h. But it is a continuous map on a subset 
of a finite dimensional space and hence conpact • So, b^^ 
Schauder's fixed point theorem (Theorem 1 *3.1l), there exists an 
E such that e(E) = E, Further, 

j 

Mel I „ <5 < K-i A 4 < m < ri-l, 

L (H^*) ~ - - 

where depends on K^, K^, K^q, and K^ 2 * 

Corollary 2 * 7 * 3 * The following estimate holds for e 

( 2.7.16 ) 1 I el 1 „ . < K. ^ h"^^, 4 < m < r+l, 

L (H-^) ~ ~ ^ 

where K^g d':pGnds on and K^. 

Proof. Replacing E by e in (2.7.14) and using estimate (2.7.13), 


we obtain 
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(2*7-17) 


iiy 


L (L^) 


2 , 1 ^15 


m—l 


Further by the inequality (1.3*14), 

I Ut I < 1 1 ^1 I, 
and from (2-7*17), we have 


iW 


L~(h1) ' 15 

Besides^we get from (2*6*25) 


< K(K_) h"*”!* 


L (H^) 


Is ^ 


1771 I 


l“(h1) 


+ lUll 


l“"(h1) 


< Kg h"'*"! I lui 1 


L (H^) 15 


Therefore, the result follows* 


We are now looking for a pair of approximate solutions 
of {U,S}, where {U, S} is the solution of the original problem 
(2 -2*1) - (2*2*4)* The Galerkin approximations and are 
g iven by 


( 2*7*18 ) 
where 

( 2*7*19 ) 


U^( y , t ) = u^( X, t ) ■» 

Sj^( t ) SS Sy^( t) f 

y = Sj^(t)x 


and Sj^, are given by the equations (2*4*4), (2*4*5) 
respectively* 
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The conclusive theorem concerning the error estimates in 
U~U f T is given as follows : 

T heorem 2«7«4* ' Under the assumptions of Theorem 2»7«3/ ( 2* 2*5 ) 

iN> ro 

condition and R^/ the following error estimates hold for 
r > 3 

(2.7.20) 113-3,11 ^ zs O(h^), 

L (0,T^) 

( 2 . 7 . 21 ) I IT-T, 1 1 ^ = O(h^) 

^ L (0,T^) 

and 


(2.7.22) 1 1 lu-u^m ^ ^ = o(h^^~^), j =1,2, 

L (0,T^/hJ(q(t))) 


whore the last norm is understood in the following senso 


( 2 . 7 . 23 ) lll^Plli 


00 = sup { I 1 <P1 1 . ^ } , 

L (O.T fH-'(Q(T) ) ) 0<T<T^ H-'(Q('r)) 


with 


( 2.7.24) 


q(t) -= ( 0 , min (S(t),S,(t))) 


Proof . Lot e^^ = ( 2*3.11 ) and (2.4.4), we get 


de. 


( 2 . 7 . 25 ) * — a(o) [e^(l,t)s + u^(l)e 2 ^j 


Integrating (2*7*25) with respect to t, we have 

t' T- 

le^(t^)l s! l-a(o) / [ e^(l,t)s(t )+u^(l)e^(t) ] dti 
t ^ 

< ; (le^(l,t)l !s(t)l + lu^(l) 1 lej^(t)l)dt 
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Since 


lujj(l,t)l < I lu^ II 

X ' — XX 


i "=xx" + l'"xx' 


< < 2 ^ 2 # for sufficiently small h and 


r > 3, 


we apply Gron wallas Lemma (Lemma I* 3 *9) to aet 


(2*7.26) Me, II 


. II „ < K(IC,K )(11T) (1,.) II 

L (0,T) ^ ^ ^ L^(0,T) 

+ I I t^(l, .)! I ) GXp (2K-K_T) 

^ l2(o,t) 2 1 


But, from (1.3.18) we have 


(2.7.27) IiL,(l,*.)ll , < >> 1 o o 

^ L^'(0,T) L^(0,T#L"‘) 


For the estimate 1 I t 


2 / choose s in the inequality 


- VJ. II >9 O, * 

L^(0,TyL'^) 

( 2.7.12) so that 2a~l = 40S and r<^lacing E by e, we get by 
( 2.6.14), (2.6.25), (2.6.32) and (2.7.13) 


(2.7.23) I 1 ' 9 9v * 

l2(o,T/L"^) 

From ( 2.7 * 26) - (2 *7.23) and (2.6.26), we obtain as estimate 

(2.7*29) Me^ll^ ^O(h^), r>3* 

L (0,T) 

The estimate (2*7-20) follows, if we note 


I IS-S, I I ^ = i Is-s, I I „ 

^ L (0,T^) ^ L (0,T) 
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Sutotracting (2»4»5) from (2-3-12) and then integrating with 
respect to t, we have 

x(t) - ^ (s^-s^) dt' - 

o 

Now from (2-7*29), the estimate (2-7*2l) can be easily derived* 
Finally, wa have 


1 1 |U~U^I I I „ 

L (o,T^; 

Hence, using (2-7-16) 
j =1/2 follows* 


H^(Q(T))) 
and ( 2 * 7 - 13 ) 


J Iu-u“l 1 


’(0/T/H^ ) 


the estimate (2-7*22), 


for 



CHAPTER 3 


NONLINEAR STEPAN PROBLEM WITH DIRICHLET BOUNDARY CONDITION 
3*1 Int roduct ion » 

In this chapter, we consider a single phase nonlinear 
Stefan problem with Dirichlet boundary conditions* Based on 
straightening the free boundary by a co-ordinate transformation, 
an H — Galerkin procedure is analysed both for continuous as 
well as discrete time Galerkin approximations* The global 
existence of the Galerkin approximation is established and a 
priori estimates are obtained* 

The organization of the chapter is as follows s In 
Section 2, the basic problem and its transformed system are 
discussed* The Section 3 deals with the weak formulation and 
Galerkin procedures* In Section 4, an auxiliary projection 
is defined and its error estimates are derived* Assuming 
existence and uniqueness of the Galerkin approximations, a 
priori error estimates for the continuous time Galerkin 
approximations arc obtained in Section 5* In Section 6, the 
global existence and uniquenvoss of the approximate solution 
are established* Finally in Section 1 , a priori error estimates 
for discrete time Galerkin ^proximations are derived* 

3*2 Problem and Co-ordinate Transformation * 

We consider the following single phase nonlinear 
Stefan problem in one spac© dimension s 
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Find a pair {U^S} so that U = U(y, T) satisfies the nonlinear 
parabolic equation 


(3*2«l) — (a(U)Uy)y = 0 in Q X (0,T^] , 

where q(t) = {y ; o < 3 KS(t)} and 0 < t with the initial 

and boTondary conditions 


( 3 * 2 * 2 ) UCy^O) s g(y)» y e l ; 

, , U(0,^) = 0 , 0<T<T/ 

( 3 « 2 . 3 ) “ ° 

U(sCT),T) =0, 0 < T < 

where the free boundaiiry S = S(t) satisfies an additional condition 
of the form ; 


(3*2*4) S-j + a(u)U 1 =0, 0 < T < T 

^'y=S(T) 

with S(o) = 1* 

Here, we malce several assumptions, which will be collectively 
called Condition B 

( 3 * 2 * 5 ) CONDITION B» 

( i) The pair {U, S} is a unique solution to ( 3*2«1 )~( 3*2*4) • 

(ii) For all p e R, 0 < a < a(p), where a is a positive constant* 

(iii) The function a(*), which d^ends on U only belongs to 

C (r) with continuous and bounded derivatives* Moreover, 
there is a constant > 0 such that 

lal, la^l, ia^pl < K^* 
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( iv ) The Initial funcTtion g is sufficiently smooth and 
g(o) = g(l) s 0. 

We organize our regularity hypothesis on the solution 
{U,S} of (3*2»l) ~ (3*2»4) according to the results in viiich 
they are used# We assxame 

u e l“(o,t^/h'^+ 1 (q(t))) n l'"(0,t^;h2(q(t))) 

nW*-' T))) n ^^^’'“(o^t^/h^(q(t) )), 

s e w^'“(o,T^) , 

R^* u e n W^'^(0 ,T^^H^(q(t))) n W^' T)) 

S e n ^(0,T^) r 

CO 

L et denote a bound for U and S in all of the spaces described 
in R^ and R^* 

Below, we fix the free boundary using Landau type 
t ransf ormation 

( 3-2*6) X = S^'^C T)y. 

Further# in order to decouple the resulting system, wo introduce 
an additional transformation in time seals given by 

(3*2*7) t = t(T) =/ dt' * 

o 

An easy calculation shows as in 2*3 that the function 
u(x,t) = U(y,t) satisfies 
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(3*2*8) ii^“(a(u)u^)^ = -a(u(l ))u^(l )xu^ ( x,t) six (0 ,t] ^ 

with the initial and boundary conditions, 

( 3 * 2 * 9 ) u(x, o) = g(x) / xs I; 

(3*2*10) u(0,t) = u(l,t) =0, t > 0 

and the function s(t) = s(i) satisfies 


( 3 * 2 * 11 ) St- + a(u(l)) u^(l)s = 0, t > O 

with s ( 0 ) = 1 • 

Here, t = T corresponds to ”*• = Note that all the smoothness 

assumptions and R 2 on the pair £U, S} carry over to {u,s} 
with the bound K 2 ( say) * For convenience, let us call these 
regularity conditions on {u, s} R^ and R 2 * ■^s u(l) =0 

a(u(l)) = a(0), a Icnown constant independent of u* 

Further, the integral (3*2*7) can be rewritten as 


( 3 . 2 * 12 ) g = s^(t), t>0 

with "^(O) s= 0« 


3*3 Weak Formulation and Gal erkin, Procedures, * 

Let H^'(I) n = £u e H^(I) : u(0) = u(l) = 0}. 
Multiplying both the sides of (3*2*3) by '"''^xx integrating 

by parts the first term with respect to x# we get 


( 3 , 3 * 1 ) + 


xx)=^(0)^x(l)(x^x^^xx 


),V e H^(I) 
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If u e l“(0/T/H^(i) n hJ(i)) satisfies (3*3*l) and initial 
condition (3«2*9), then it is called a weak solution of (3«2*8) - 
(3*2-lO). 

H^>^Galerkin procedure . Let ^ c fl be a finite dimensional 
subspace^ satisfying the approximation property (l*3«i6), for 
k = 2 and inverse property (l»3*'9)* Then^ the Galerkin 
procedure is defined as follows : Find u^ : (0 ,t] -*• such 
that 

(3*3.2) (n^3^X5^) + ((a(u^)uJj)^,X^)=a(0)uJj(l)(xu^/X^^), 

•u 

and u (X/O) = Qj^g(x)/ 

where Qj^ is an appropriate projection to be defined later* 
Moreover/ the Galerkin approximations Sj^ and of s and x 
respectively are defined by 

ds. V, 

(3.3.3) ^ = *-a(0) u^(l) Sj^, t > O 

with Sj^(O) =1 

and 

dx 

( 3 . 3 . 4 ) - Sj^(t), t > O 

with ~ 

Next/ we define a fully discrete, 0(At) - correct 
method for approximating {U/S} based on backward differencing 
in time. Let At > 0, N = T/At eX and t^ « nAt/ n e ^ . Also 
let = <p(x/t^)/ and d^. = (fP )/At* Denote the 
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approximation of s by W and u by Z. Assuming and z” to be 
known, we determine and as follows : 

( 3 .3 .5 ) (<at25,X^)+( ( a( z") )^,x^)=a(0 )zg( 1 ) ( 

with Z*^ defined approximately 
and 

(3-3*6) = -a(0)2^*‘'^(l) w"*'’! 

with W° = 1 * 

Moreover, the discrete Galerkin form of (3-2#12) is given by 
( 3-3-7) "^h “ (W""^^)^’ 

with = 0- 

3*4 Associated Projection and Related Estimates - 
Set 

(3-4-1) A(u/v,w)=( (a(u) Vj^)^,w^j^)-a(0)u^( 1) (xv^w^^),for 

u e W'-'^’d), V and w e H^(l)«, 

For the operator A the following Lemma holds » 

Lemma 3-4-1 (i) There is a constant M > 0, depending on 

K and I lu 1 1 such that 

^ L 

( 3-4#2) lA(u/v,w)l < Ml Ivl 1 ^Mwl 1 2 # u s W^‘'*”/v and w e 

(ii) Garding type inequality holds that is; there 
exist constants ot and P such that 
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(3i4ii3) A(u/viv) > a I Ivl I 2 -* PI Ivl I for u B and 

2 1, 

V 6 H (l) n H^(l), where ^ may depend 5n and ) Ju^l I ^ 

L 

Proof I! From (3»'44l)^ we get 


|A(u/v,-w)l < Ha(u)v^^iW^^)l + l(aj^(u)u^v^w^^) l+la(0)u^(l) 

< Viiv^ii Iiv^ii+nvi^jiv^n nw^ii+iu^(i)i iiv^i i 


' '"xx' ' 


< J llvlUHwllj* 

L 

For ( ii) t we have 


A(u#v*v) > (a(u)v^^v^^)-{ t+la(o)u^Cl)(xv^,v^)l} 


> ai l lUj^'i I „ 

Xj 

Now applying the inequality ah < 
term and noting N^'ll < 


I )v^l I 



for '<p e 


'''xx" 




2 

Eb for the second 

we get the requircsd 


result* 


Let u e Sj^ be the projection of u with respect to the 
form A that is 


( 3*4*4^ A( X ) is 0^ X 6 

We shall consider the existence and uniqueness of u subsequently. 

Let l(u) and L*(u) be the one dimensional elliptic 
Operator and formal adjoint associated with (3'*4*4)jf where 
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(3*4*5) L(u)w = 1^ (a(u)vrj^) - a(u(l ))u^(l) 
and 

(3*4*6) L''''(u)w = I— (a(u)w) + a(u(l)) u„(l)(xw)„* 

Standard regularity theory combined with uniqueness results, 

see Douglas et al* [ 13 ] modified for one dimensional case 

imply that the Dirichlet problem for L has atmost one weak 
•1 

solution in H‘^(l)* Following Douglas [ 16 ] v- one shows that for 
5! s L'^(l) there exists a unique solution of the Dirichlet problem 

L'“*(u)g = 4 , X s I 

g a 0, X = 0 and x = l|i 
satisfying the following regularity condition 


I I g 1 I 2 < i I W I I / 

where the constant depends on u and its derivatives* 
For 'ie L^-(I), define cp e H^( l) D hJ( I) by 

( 3*4*7) L*(u) ^ ® ^ 




0 , 


By regularity assumption and boundary condition (3*4*7) (ii) 


Since s H^( I)* the regularity of the Dirichlet problem for 
o 

the Laplace operator implies that 
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(3*4*8) n cpI 1 ^ < C^l I ^1 1 * 

Integration by parts shows that 

(3*4-9) (v^,L'"(u)<P^^) = (L(u)v^^^, (p^^) 


= A(u;v 3 ^,<p),for e H^Cl) fi H^d)* 

Now, with the above preparation, we show the existence of u 
defined in (3-4-4) , see Douglas [l$]- 

o » 

Lemma 3*4-2- There exists a unique solution u s of (3*4-4,/, 
atleast for sufficiently small h- 

Proof - It is sufficient to show uniqueness- Let u^^ and u^ 
two solution of (3-4-4) and set 0 = ‘u^— u^# Let <P S fl 
satisfy 


L'*Hu) <P 


XX 


-^XX- ^ ^ 


V. 


XX 


x=0“ 



0. 


Multiplying both the sides by 0 and integrating by parts, we 
use ( 3 - 4 * 9 ) to get 


110x11 


Since A(u 70, X) 
we obtain 


2 = (e,L»(u)=p^) = (i.(u)e^^ 

= A(u;0,<P) - 

= 0, X e and 110 11 < 1 10^1 1, for 0 e 


I 101 1^ < A(uf0,<P-X) 


< MI 101 1 2 


inf 1 I<p-Xl 1 , 

Xegj, 
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< K^Mh^i lei i 2I l<i>l 14* 

Applying regularity property that is 1 1 cpl 1 4 4 j 10 I I 2 we get 

1 10 1 1^ < K^c^Mh^i lei I 2 * 

Since 0 8 A(u/0,0) ss 0 and by Carding type inequality (3*4*3) 


Thus* 


al 101 1 ^ < PI 101 I 2 

< K^C^M Ph^l 10 1 l2* 


0 = 0 * for sufficiently small h. 


Hence the result* 

cv 

Let T) = u*^u» At first we obtain some estimates of V 
and its tenporal derivative Tl^ for our later use# 

The following Lerrma will prove very convenient for that purpose# 

An analogous result was used by Arnold and Douglas [ 4 ] in the 
2 

context of L -Galerkin procedure# 

Lemma 3#4#3* Let 5 6 H^( l) fl f^(l) and satisfy 

(3#4.10) A(uf?*X) = F(X)* X e 

where F : H^( l) fl H^lC l) -+ R and linear* Let there exist 

o 

constants and M 2 with M^ > M 2 , such that 

IP(<P)| < M^^l icpl I 2 # (P e H^(I) n H^C I) 

( 3 #4 # 11 ) 3 .nd. 

1F((P)1 < 1<PI I 4 / <P e 
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whore h'*(i) = Ju c H*(l) n I^CDiU^I =u^| =0}. Then, 


for sufficientay small h 


ll$ll j < K [m + inf II§-xll J 

^ XeSj^ ^ 


( 3*4*12) and 


1 1 §1 I < [ (M^+ inf I 1 $-Xl I 2 )h\M 2 ] t 

where = K 3 (ct, P, is used as a generic constant* 

Proof * First note that 


A(u?$,f) = A(uj§/§-X) - F($-.x) + F(§) 


< Ml I i 1 I 2 inf I I $-X I 1 2+M. inf f I ^X1 1 +M I i §1 I 

o 

By Carding type inequality ( 3*4*3) , we have 

CtIl§ll2~PlI§lI^< (M!1§11 2+M ) inf I I Xl 1 2 + ^3. ^ ® * * 2 * 

XeSh 

Using the interpolation inequality for that is 

1 1 § I 1 5 , Mill I I § I I 2/ we get 

( 3 * 4 * 13 ) 1 1$!!,, < (a)~^ [pll#|l+M inf 1 r#-Xl 1 „+(K +l)M,]. 

Z Q 2 0 1-^ 

2 

For L -'estimates, let us recall Aibin— Nitsche^ s duality argument* 
For a given 'i e L^(l), define cp e H^( l) by ( 3*4*7) * Multiplying 
both the sides by §, we obtain 

(#,15;) = (#, = A(u/#, (p) 
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= inf [a(u/§,<P-X) F{q)-,X) + F(q))] 

XeSh 

. < (Ml |$1 I 2 +M^) Inf 1 l<P-X! 1 2 + M 2 I 

Xes. 

h 

< [(Ml I §1 1 2 -!-Mj_)KQh^ + M 2 ] I 1 <Pl 
Applying regularity property (3<»4#8)j^ we have 

($,^') < [k^(M 1 Ifl l2+M^)h^-fM2] ll«IU 

Thus^. 

(3.4,14) 1 I $11 < [Kq(MI I$1 I 2 + M^)h^ + M2]^ 

Substituting (3i^4#14) in (3*4 . 13 ), we get 

II $11 2 < (a)"^ [P{C^Kq(M1 1 $1 1 2 +M^)h\c^M 2 }+M inf 1 i $- Xl I 2 

Xeih 

t (k^+i)mJ . 

For sufficiently small h, we get the required estimate 
(3»4*>12) ( i) for Il$ll 2 * 

Now^ making use of 1l$ll2 in (3*4*14)< we obtain the estimate 
for 1 I $1 I , 

Remark* Since $ e fl the estimate II $11. can be found 
out either by interpolation inequality between I I$1 I and 1 1$| I 2 
or directly substituting $^^ for 'i and using regularity 
I i(pi I 4 < c^i i$n 2 ’ 

The next Theorem contains the error analysis of n and for 
the associated projection (3*4*4)* 
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Thfsorem 3»4*4*- For t s the error V = u— u satisfies 


( 3*4-15) 1 iTIllj. < lulij^, j = 0,1,2, 2 _ m < i+l 

and 

(3-4-16) 1 IHt* I j < I Ijn+I lul Ij^), j =0,1,2, 2< m < r+l. 


where the constants K* = K,(K .K_) and K- = K_(K ,K. ,K_,K., 

^ 4 O' 3 5 5 O' 3 4^ 




Proof* To get the estimate (3-4-15), for j = 0 and 2, put 
i = T) and P s 0 in the previous Lemma 3*4*3, then 


I 17 ? 11 ^<K 3 inf Iir?~xil 2 
xsih 

o 

< K 3 inf I lu~vl I 2 / for v = X + u e 

VSSh 

< h"'”^ 1 lul , for 2 < m < rtl* 

Similarly, 

1 fT)l 1 < 2 < m < rhl* 

For j = 1/ the result (3-4-15) follows by interpolation- 
To estimate T)^, we differentiate (3-4-4) with respect to 't' 
and obtain, 

(3(.4*17) A(u/h^,X) = "(l^j ( [^ a(u)] t)^) , X^) 

+ a( 0 ) u^(l) (x ri^ Xxx? * 

Identifying the right hand side of (3 ’4*1 7 ) with F(X)/ we note 
that for ^ e n 
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(3.4,18) IP(<P)1 < Kg) inll^ 1 l<Pll2/ 


where Kg depends on K , 1 lu 1 1 ^ , I |u 1 I ^ , 

L (L ) ^ L (L ) 

O 4 

Further, when <P e H (l), integrating by pairts the right hand 
side of (3'»4*17) for X = tp we obtain 


F«p) = (a(u)t ■n^ 

= -(>l,(a(u)^ ^xxx’x' - iP^xy ' 

So, 

( 3 . 4 . 19 ) 1 F(CP)| < Kgl IT?! I ||(pl)^. 

Thus, Lemma 3.4»3 is applicable and we get 

1 inj.! I2 < K3 [Kgl l77l I2 + inf 1 I'D^.-Xl I2 ] 

XsS, 

Here, \<re have used the estimate (3.4*15) for j = 2 and the 
approximation property 

inf 1 l |5 «• XI l2 = inf • I2 < 2 < m < r+l. 

XeSh eSj^ 

Further, using Lemma 3*4.3 and the estimate (3»4*15) for j = 0 
we get 

1 II < [ ( Kg I I T? I I 2+ inf I I Tl^- Xl i 2 + Kg I I pi I ] , 

‘xe§h 


Hence, the required estimate for llp^ll is obtained# 
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For j • = If the estimate follows as usual by interpolation* 

We shall also need later the following estimate for 

T}j5,(X,t) . 

Th^iorem 3*4*5* There exists a constant = K^(a/M,K^,K^) 
such that 

(3-4*20) It)j^(l,t)l < I lul 1^/ 2<m<rtl. 

Proof* Let us define an auxiliary function <P 6 H (l) fl H^( l) 
as a solution of the following problem 



Multiplyino both the sides of the first equation by tj and 
integrating by parts# we get 


laCo) -n^Cl/t)! = lA(uyn, (p)l 

= 1A(u/7?,(p*‘X^ J 


< Ml 1 7)1 1 „ inf 1 l<p-xl I o 

^ o * 

XBS^ 

< K^Mh™''^'! Id H 2 ^ • m* 

Now applying the estimate (3*4*15) for j = 2 and noting 
•a( •) > 0.t we have 
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irf^(l»t) 1 


a'^^K^MK.h 
O 4 


2 (m- 2 ) 


l<»>l 




2 <, m < r+l/ 


which is the required result* 

3 ‘*5 A Priori Error .Estiroat^s for Continuous ^iroe Galorkin 
AqT>r.dximatioh s » " 

Throughout this section we assume that there are 

, h o 

constants K'' and h^ such that a Galerkin approximation u e 

defined by (3)«3*2) exists and satisfies 


( 3*5.-l) I 'lu^i 1^0 < for 0 < h ■< h . 

L (H^'(P)) - ° 

It. 

\vhere u (x^O) = Qj^g 

o 

and A(g;g~Qj^g.^X) =s 0^ X e S^* 

Clearly.> U^(x>0) S u(x^O)» 


We use below the standard notations-* 

It li 

e = u~u - u '^u and 71 s= u— u.^ then e « T) & * 

‘ro h 

.' t.heorem 3*5'vl* Suppose ‘H t= u-u satisfies '( 3-*4'*4'i) and u /, 
defined by C 3’*3!»-2) ' is the Galerkin approximation of u^* Assume 
further that ( 3‘>5'**1 ') hOldsj/ then there is a -constant 
Kg = K^(ar,A/K'/,I<;j^vK^,,K^/,K. 7 )su'ch that 


(3.5'.-2) ^ .'tal 

L 




-'+ ^Muil 3 < m 'itl-* 


.Proof* From ( -3"’*4-"*’’4) and (3>3'*'l) with 'v = .‘X /<* we get 
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Subtract it from (3*3*2) to obtain 

( 3 * 5 * 3 ) (^.^/>^) + A(u^?u^,X) - A(u/u,X) = * 

But 

A(u^fu^, X)''A(ufu, x) =A(u/ if X a(u)"-a(u^) )u^ )^^) 

~ ( (a(u)-a(u^))u^^ ^xx^* 

Substituting this in (3*5*3) and setting x = we have 
(3*5.4) 11 ^ IU^Il^+A(un,S)=(’ltx"^,)+(| 5 j(a(u)~a(u^))u^ 1 ^) 

+ ( (a(u)~a(u^))ujj^, l^^)-a(o) ^xx^ 

+ a(0) 

Applying Carding type inequality (3*4i»3) for A(u;S/S) and (3*5.1)/ 
we integrate by parts the first term in the right hand side of 
( 3*5 *4) to get 

I ^1 1 ^11 \al U1 1 2~ PI i ij < 1 ’ W K2)(limi^+ 

\lil\^)\\i^i\-^K^K‘H\{l)\ 1 U^^N+K^K'^1 yi)! 

V y O 

Using ( 1 * 3 * 15 ) for lx(l) and I Ixl 1 < IIX^^ll/ for Xs S^/ we 
have 

(3*5*5) ^il 2lI^+2Sl 1 21 1 ^2P lull 1+21 )T)^1 I lUn2+2Kg( I l7?i l^+r]^!)) 
Will 2+4Kg 1 Ul li N I+K^K*! 1 21 ly ^1 1 il I 
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Prom which,we obtain by using ab < 4.^-1 

a,b,e > o) 

3=E lISlll + &\\ll\l < Kj^o(IlT)tli^+ll’7M^+|T?^(i)l2)+ic^^(6)ll gn2 

+ (Kj^q + 2P) 11 l\\^. 

Now with appropriate choice of 6, 

fVJ 

to g* With this choice of e, we get by integrating with 
respect to t^ from 0 to t and using Gronwall's Lemma (Lemma 
1*3*9) the following 


1 1 SI l^(t)+a / 1 1 SI I^dt' < exp [(K^Q+2P)t] / ( j iTi^J I |pl 1 


+ in^(l)l^') dt' , 

where depends on and K 2 * 

Now taking supremum over all t e [o,t] and using the estimates 
( 3 * 4 * 15 )/ ( 3 •4*16) and (3-4 *20), we obtain 


SI 


. (H^) l2(h 2)~ ^0 5 t l2(h™) l2(H^) 


K^h 


m~l , 


lul 




2(m-2) 


1 lul I 


L^(H^) 


Therefore, for 2(nv~2) > m~l that is for m > 3, we get the 
desired estimate (3*5*2)* 


Corollary 3*5 *2* Let the inverse hypothesis ( 1*3*9) be 

o 

satisfied for X e alongwith the assumptions in the previous 
Theorem* Then, 
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(3*5*6) lUII^ ^ <K_h"^^[j]u. 


+nuii 


L (H^) ^ L2(h"^) 


where the constant - K^^Ck^^Kq) 


KJ 

Proof * For ^ e we get by inverse hypothesis (l»3*9) 


' * W ‘ 00, 2^ 1 CO o • 

L (L^) ° ^ L (L^) 

Thus, (3«5*6) follows immediately using the estimates (1*3*14) 
and ( 3*5 • 2 ) » 


Prom Thcsorems 3*4*4, 3*5*1 and Corollary 3*5*2, we 
get the following Theorem* 


Theorem 3*5 *3* Let the solution u of (3*2*8) ~ (3*2*10) be 
sufficiently amooth to satisfy the regularity conditions 
Further, suppose that there are positive constants and 
K*(K* > 2K) such that an approximate solution u^ e ^ satisfying 
( 3 * 5 * 1 ) exists in I X (0 ,t], for 0 < h < h^* Then the following 
estimates hold for r > 2 


( 3*5*7) 
and 


llel 


’(H^) “ 


K. 


12 


(3.5*8) Hell 5 < K , 

L (H^) 

Where K^2 = ^12^ ^2' ^4'^^ and ^4'^11 ^ constants. 

Further, for sufficiently small h and r > 2 


( 3*5 *9 ) 1 lu^l 

and consequently, 


'00, 2 ^ 2K2 < 

L (Hd 


K ^2 es well as 


can be chosen ind^ondent 


of K*. 
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Proof* The estimates (3*5«7) and (s+S^S) are iirtnediate from 
Theorems 3 *4 *4, 3*5 •! and CJorollaxy 3*5 *2 by triangle 
in equality* 

To see (3 *5 •9), we note 


I u 


* CO . O . 

h (H‘‘) 


< I lul I 




+ I lell 


“(h2) 


< + K^3 


< '^ 2 * sufficiently small h and r > 2* 


Finally^ the Galerkin approximation of the solution 
U(y, of the original problem can be defined as 

(3. 5 •10) U^(y,T) = u^(x,t) 

where 

( 3.5*ll) y = X 

( 3-5-12) T = 

where Sj^ and are defined by (3*3*3) and (3*3*4) respectively. 

Theorem 3*5*4* Under the assumptions of Theorem 3*5 *3^ (3*2*5) 
condition B and the regularity hypothesis R^, the following 
estimates hold for r > 2 


( 3 *5 *1 3 ) 


1 1 S-Sv, 1 i „ 


= o(h^); 


1 IT - Thl 1 oo , = 

^ L ((0,T^)) 


( 3 *,5*14) 
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and 


( 3 * 5 * 15 ) IIIU-U^MJ 


l”(0,T /H^( IXt))) 


= j =1,2. 


where the last norm is understood as in (2*7*23) * (2*7«24)'. 


P roof * Subtracting (3*3*3) from -(3*2*11) and then integrating 
with respect to t, we get 

t t . 

is-sj^i < K, [;( in (i)i+u (i)i)isidt'-(- ; iu^(i)i is-s>,i dt*]. 
o o 

By applications of Gronwall's Lemma (Lemma 1*3*9) and the 
inequality (l*3»16) for §^(l)/ we get using the estimates 
( 3 * 4 * 20 ) for m = rfl and (3*5*2), 


< h , for r > 2* 

The estimate (3*5*13) follows immediately, if we note that 


‘ >S~Sj^|] « IIs-Sj^ll • 

" L (0,T^) ^ L (0,T) 

f 

Moreover, the estimate (3*5*14) follows easily using (3*5*13) 


Finally, we have 
I I lU-U^ l 


l“(o,T ;H^(Q('i:))) 


< I lu-u f I 


l“’(0,T/H^(I)) 


Hence., using Theorem 3*5*3ir the required estimate for j =1,2 


follows * 
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3*6 G lobal Existence and Unlcrueness of tho Galorkln 
Approximation « 

Nov7 wg consider the problem of eocistonce of the 
Galorkin approximation u in the domain of existence of u-. 
Towards this end, let us recall (3*5*3), (3*5*4) and got 


( 3*6.1) ( 

+ ([a(u)-a(u^)] u^^,X^x^~a(0) ^^^l) ( xu^/X^) 

+ a(o) ^^( l) ( 3cu^,X^^) • 

But 

•I 

/ % / hs , 3a(u~gG) ^ jt 

a(u)-a(u ) = cd^. 

Replacing u^ by u~e in (3*6.1) and with minor adjustments, 
we obtain 


( 3*6.2) ( ^^/X2^)+A( u/ ^^xx~^xx^ ^x^ 




o 3u‘ 


(u^-'e^)d^ V ^ ^^^^'^x~°x^'^xx^ 


'XX 


+ a(0)yi)(x(u^'-c^),X^^). 


Substituting e by E(x,t) for some function E e it 

follows that 

(3*6.3) ( 
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XX 




E^)dS, 


XX 


1 

)-(; 
o 


3^a(u-6E) 

3U^ 


+ a(o)yi)Cx(u^-E^),Xxx>/ X ^ V 

This is a lincrr ordinary different ial equation in Therefore^ 
for any E = E(x, t) there exists a unique solution i of (3«6#3) 
with • 


( 3*6*4) ?(x,0) = 0 

in the interval (O^t], 

The equation ( 3 * 6 * 3 ) defines an operator such that I = *^E 
2 1 

for each E e H fl Since e = n-?# therefore 

( 3 . 6 . 5 ) e ='n-?7E^ for each E e 0 hJ * 

To show the existence of a solution u for (3.3*2)/ we need to 
show that the operator equation ( 3 * 6 * 5 ) has a fixed point. 

In otherwords, we are looking for an e( e) such that e(E) = E* 

Theorem 3.6*1* Assume that the finite element space 

satisfies the inverse property (1*3*9) and u is the unique 

solution of ( 3 * 2 * 8 ) - ( 3 . 2 . 10 )* More over,let the condition 

B ( 3 * 2 * 5 ) and the regularity hypothesis for u be satisfied* 

Then for 6 > 0, there exists a solution u e Sj^ of (3*3.2) 

satisfying 1 lu^u I 1 «> o, 5. <5*' 

L (0,T/H^) 
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Choose X =s: 2 in (3 #6 # 3 ) to ottain 

I Ie Il^n\aii?:ii2«piig||2 < II^^I, lis^li+ic^ciinil ^+l\l\\ J 

L Lf 

(llu^j^ll+IIE^II)l|J^^l|+K^(||„^|| + ||!^||)C||Uj^|l „+nE^II J 

L L 

1 1 ^3,^1 i+K^( iiu^i I ^+i\E^\\ J^( I mi i+ii Dn^^i i+K in^{i)i 

L L 

I 

( I lu^l 1 + 1 lE^I-l )I1 ^^^1 |+K^( I fu^l I+! lE^l I )u^(l) I II I . 

Using (l*3*'lo) for and applying Young' s inequality for 

this term and some other terms, we get 

§:£ I I SI lj+2al Ul i2 < K^gCe)!! SI 1 2+Kj^7(s,K^,K2){ 1 in^n^+ 

( 1 Ihl I Ihl l^)(l+l !E| l2)+Mhl 1 ^(1+21 I El l^Ci+l lEl 1^) ) + 

L 

1 \(l) I ^ (Itl lEl l^)}+K^Q(e,K^,K2)£ I I El 1 ^(l+I lEI Ip+P] I I SI 

Choosing e appropriately so that = 2a, integrating with 

respect to t and then apnlying Gronwall's Leirma (Lemma 1#3^9), we 
obtain 

! ISIl^Ct) < ejq> [K^g{llEl|^„ ^ il-i-] ] E\ ^ t“^t’ ‘ ^ 

+ imiL^-(l+llEll^„ „ )+|n (l)l^(l+l lEl 1 ^^ ) 

^ L (H"^) L (H"^) 

+ llhl!^{l+2l lEll^^ .(l+llEll^^ 2 )} ] dt'. 

L (H^) L (H^) 
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From Theorem. 3*4 *4 and Theorem 3*4*5, we get 
(3.6*6) lUIl ^ ^^19 [h^'*'^+(l+l I Ell _ 


L°*(H^) ~ 


l“(h2) 


+ h^+^{l+l|Ell ^ (1+llElI ^ )}] , 

L (h2) L~(h2) ^ 

where K^g depends on K^g, 1 1 E| I ^ 2 ^ * 

L ( H ) 


Thus, we have. 


(3.6*7) Hell eo 0 < ll^" CO 0 + 11^" oo 9 , 

h iU^) "■ L~(H^') L (H^) 


-li 


< 11T71I CO 2s l^' 

L (H^) ° L (H^) 


( 

Now, for I I El I „ /j <6 we get from (3.6*6) and (3*6*7) 
L (H^) 

"ell „ . < K 0 , 

L (H^*) ^ 

where K^q = 

Therefore, for sufficiently small h and r > 2 


It is easily verified, see Theorem 2.7.2 that the operator 
is continuous on finite dimensional space and hence compact* 

So the operator E -* e has a fixed point by Schauder's fixed 
point theorem (Theorem 1.3*11). Hence the result. 

The uniqueness can be proved following the argiaments 
in Theorem 2*5*2, except that the estimate of <}) 2 ^(l,t) is 
given by the (i.3.l6) instead of (1*3*15) » 
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We formalize the above in the form of a Theorem* 

Th eorem 3*6*2# Let u satisfy the regularity condition and 

h o 

let K > 0* Then there exists one and only one solution u e 

of (3*3*2) in the ball { I lu-u^l i - < for sufficiently 

L (H'^) 

small h and r > 2* 


3*7 A Priori Error Estimates for the Discrete Time Galerkln 
Approximations* tt,.,. ^ - 

In this Section, we first develop a priori bounds for 
the error u^~Z^ and s^-W^ • The results obtained will be similar 
to those obtained for the continuous case in Section 3*5 * We 
assume that condition B (3*2*5) and the regularity hypothesis 
for {u,s} hold* Let ^ = z’^-u” and e’^ == = 

Theorem 3*7*1* Let Z° be determined by 
(3*7*1) A(g;g-Z°,X) =0* 

There exist constants = K2j^(p,a,M,K^,K^/K2i-K^^Kg and K^) 
and h^ > 0 such that if h < h^, t = o(h) and r > 2 

( 3.7*2) sup I U1 I ?+3At I I I I 2 < ^ 21 ^ (At • 

^n n=0 

t 

Proof * For t = t’^'*'^, the projection (3*4*4) is given by 
(3*7*3) A(u'^’^^/’f‘^^,X) =0, X e Sj^* 

From (3*7*3) and (3-3*1) with v = X , for t = t^"’’^ we get 
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3u 


(3.7.4) (atVXx>+A(u"'^^*u’’'^’',)0= -(at 


Subtracting (3»7.4) from (3-3*5), we obtain 

(3.7.5) (dti“,Xjj)+A(z"72”+^,X)-A(u‘’'^^/u”'^'',X) 

“ (at n^x,) - (at“x-S-^)'^^>«x)- 

But, 

(3.7.6) A(z"/z'^+^,X)-A(u”+^*u’"‘*'^,X>=A(u''+^/^”+^X)-.(|:;;: ( [a(u'^)-^ 

gx • 

a(z’^)] Zx'^^)/Xxx^'-a(0)eJ(l)(xz5^’^^,Xx.3^)-(|^( a(u“+^) 

-k<u’^)]z^‘^^),X^)-a(0) At Ot • 

Integrating by parts the last two terms on the right hand side 
of (3*7.5) and using (3*7*6) in (3*7.5), we get by choosing 

X = 

(3.7.7) (dt^s5'^M+A(u”+^*s“+^f'"^)= -(at’).^^)+(at>i"-(|^)”'"^ 

5j;«)+([a(u”)-a(z‘‘)] 

t([a^(u'')-a„(z'‘)] 

-a( 0 ) 4( 1 ) ( x^l. M [a(u"«) -a(u" ) ] ) 

et’->+a(0) At At '7 x<iHx25'*'^.^^>* 
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Sinca, i I 1^) ana 

A(n>'«,£"-l,jn-HX, > a,U«l,, 2 . pH,"«n 2 ( 3 .,.,) 

(3.7.8) ^ c I I ^+1, , 2_| , J.n| I 2)^1 1 jn+l|.| 2 _p| , jiM-l, 

<lld^r,"|| lli^ln + lla^^l, ||j 5 «u+Kjll 25 «n 
( 1 1 ^"l I +1 |T)"| 1)11 ^i| l+ic ( I i ^1 l+l I rPi I )| |z”'''l|| 

L 

ll!5+’-M+K^(IU"||+|m“||)M,iJll „il^+l|I „ll^^ll 

L L 

+ K^ln;j(l)l I|Z5+1|I M^l|H.KjliJJ(l)| ||Z"1|| 
il2^«ll+K^ Atllatu'’ll iH"t^ll+K, itlla^u”!! 

J-* L L 

ii^+^ii+Ki AtiatuJ(i)r 1125 + 1)1 1155 I 11 , 

where -. (|g)^‘+l. 

Using generalized Young's inequality (Inequality III Of Chapter 1) 
forK^IJ^d)! II25+I11 I155I1I. using l^d) I < 11^1 1+11^1)1/2 

' ' ^x' “'1 ine<iuality ab < ^ + | for remaining terms 

we get 


C 3 . 7-9 ) ^ (I I I 2^1 I I 2)^^| I ^n+1, J 2 ^ pj J ^n+lj J 2 ^ I J J ^n, , 

+ (7 + |) e II ^■^^llJ+K(K3^,K2,e)nz^^^J i2(h.)J2^-!-1] |2) , 


to to 
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I 1 l^+K(Kj^,K2*e)llZ^+^) |2( 1177^1 |2+in^(i)j2)+ 

K(s) 1 Iclj. n”l 1 ^+K(K^,K 2 /e) At^i I 1 1 ^«,+ 

L 

+ K(e) llajjll^. 

L 

We next multiply each of the terms of (3 *7 #9) by At and sum on 
n, n = 0, to obtain 

(3.7.10) -I llg'^IlJ+Ata s 1 1 I? < Atp s I I l?+Atp[ 1 ^11 J 

^ ^ n=0 2 - 1 1 

+ AtK{K^,K2,e) I i l^Cl+II z’^’*‘^l I 2 ) II ij 

+ (7 + ^)e At £ 1 1 1?+ I At £ 1 I 

* n»0 ^ ^ n=0 ^ 

+ K(K^,K2,e) At £ 1 l^( I Ih^l 1^+f 71^(l)l^) 

+K(s) At £ I I ^At^CK /K ,e) E Ilz’^'^^llJ 

nssO nssO 

( I i 11^11 ^ '%^x^ 1 ^)+K(s) At £ ! I a 11 2. 

L L nwO 

But, 

(3.7.11) At J^lla„n2< (it)2 

( 3*7*12) 
and 

( 3 . 7 * 13 ) 


At i ll<it'i”lli < Cllu^ll^j . 

ns=0 ij \ri ) 

N^l *-» o 0 

At £^ ldj.u^(i)l < c ' '^xt* 
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Choose e and At in such a way that a • I5e > a- and 4-' Atp = a- 

Jm ^ ^ 

in (3*7.10), then we get using (3.7.11) - (3.7.13) 


(3-7.14) ttjl l^+2At a. S 1 !S”l I? < At '^E^£P+K(K.,K_)1 IZ”'*'^I|2 

n=0 ^ n=0 ^ ‘ 2 2 

(1+llZ^-^ll l2)j 1 1^^112^. K(K.,K_)At 

( I ln”l l?+ln^(l)l^+l let 77^1 l2)+(At)^(K.,K ) sup 

0<n<J 

9m0 tmm 


Now assume that there is a K* > 2 K 2 such that 

( 3 * 7 . 15 ) Ilz’^ll^^K*, for n = 0,1, ,-.,J. 

Then, applying (3.7.15)# (3.4.15)# j = 1# (3.4.16), j = 0 and 
( 3 - 4 . 20 ) we get 

aj Ig'^l 1?+At a. S Ilg”ll?<At s^{P+K(K. ,K-)k*^( 1 +K*^)} I I21 1 2 
^ -LI jj_o ^ n=0 ^ 


-I- K(Kj^,K2,K’M{h^^+h^^^^^ + K(K^,K2#K*)(At)^ 
t K2(At)^. 

Finally ^applying discrete Gronwall's Lerttna (Lemma 1.3-10) we 
get for r > 2 

(3-7.16) I?+At 3 S 112^11^ < K (K*) {h^^-f( At) - 

^ n =0 2 

Now we show that for sufficiently small h, K 2 ^ is independent 
of K*. Applying inverse property (1.3.9) we get foi* t = t , 
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Mz”!l2 < lle”ll2 + Ilu’^112 

< 1 !r)^l I 2 + il2’^ll2 + i 

< I + ^2 

< + K^h**^ (h^+At) + Kg* 

So, for suff icien-tly small h and At = o(h)/ 

(3-7-17) sup Ilz’^ll„< 2K^. 

0<n<J 

Therefore, the required result follows* 

Go roll ary 3-7*2# Suppose that all the assumptions of Theorem 
3#7-l hold and that the space and time discretizations satisfy 
the relation At = o(h). Further, suppose that satisfies 
the inverse property (l.3#9)# Then there is a constant 
KgjL^ such that for r > 2 

sup 1 1 ^ < K' j^th^^CAt)^ + - 

n 

This follows from the inverse property (1*3*9) and the 
estimate (3-7 *2)# 

The Theorem to follow shows the estimates of e and e^, where 

e- =: s~W at discrete time level t - t « 

1 

Theorem 3-7.3. There exists depenfling on 
and Kg such that for r > 2 

(3#7.18) sup ( llel I < KggI (^^ 

. n 
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Proof * The estimate for lle^ll^# for n « 0, l^«*vN follows 
from Theorems 3*7 *1# 3*4*4 and triangle- in equality# 

To find an estimate for consider (3«2*ll) for t = t and 
ottain 


(3#7.19) d^s’^ * d^s" - - a(0)u”'^^(l) s"^^# 

Subtracting (3*3#6) from (3#7#19) and multiplying the resulting 
n+1 

one by e^ , we get 

(3*7*20) <dj.e^,€^^> s <d^s’^-(^)”'^^/e^'^^>-<a(0)^”**^(l)e^ ^ 


with 


e° = 0# 


Here <u,v> » uv# Now the left hand side of (3*7 *20) is estimated 
by 

( 3 . 7.20 

Prom (3*7*20) - (3*7*21), we have 


where = dt-S® - 

Using Sobolev Imbedding Theorem (Proposition liS.t) for one 
dimensional domain for and applying the inequality 

ab < — + - b^, a,b > o; s > 0, we get summing on n, n=0,l, 
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(3.7*22) 


1 

2At 


J I 


J-1 « J-l , 

s Icf- r+3 E ‘'+K, E 

n-sO n=0 ^ 


J**l 

- « lz”‘^^(l)l 
•^naO ^ 


lef ll2 


j 1 

+ K(K^,Kj>) E |2 +k(K.,K„) E 

n=0 ^ 12 2 

Multiplying both the sides by At and applying the estimate 

(3.7-2) for the last term in (3.7.22)/ as well as (3.4.20) for 

T)^(l)/ we get by (3.7*17) 


( 3 . 7 . 22 ) I le'Jl^ < (At)^l Is^^l +3At E^leJl^+2K.K At 

^ l 2(0/T) n=0 1 ^2 

J-1 

E^l e^l 2+At( 3+2K,K2)1 ejl 2+K(Kj^,K2/K^)h^^ ^ 

+ K(K^/K2/K23^) h^^. 

For small At/ ~ At (3+2K^K2) = 3 ^ and now applying discrete 
Gronwall's Lenma (Lemma 1*3.10), we get the required estimate- 

The following is a direct consequence. 

Corollary 3*7.4* Suppose that all the assureptions in Theorem 
3-7*3 hold and that Sj^ satisfies the inverse property (1.3*9). 
Then there is a constant K -^3 = ^^^o'^4'^23^ such that for 
At ss o(h)/ 

sup I lell| < K'g{h’*^(At)^ + h^^2j ^ 
t^ 

We are now looking for a discrete time Galericin approximation 
of {U/S}/ where {U/S} is the solution of (3-2-1) - ( 3.2-4) - The 
approximation is given by 
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«S'yh*-'h> » V'h> = 

where xj is aeflned by ( 3 . 3 * 7 ) and - xTii”* 

Sirtiilajrly the d.igcjrefte time Galerkin approximation is 
defined by 

®h ® V''h’ = ' 

R^grna^* Here corresponds to t^ under the approximate 
transformation* Consequently «the time mesh in corresponding 
to a uniform time mesh At may not be uniform* 

Let AtP as and AT>, = max ZSr?*# 

0<n<N ^ 

Finally^ we establish the following theorem for the discrete 
time error estimates where ^ a S S(tJJ) 

and is the time cssrresponding to in the original 
t ransf or r«mit ion* 

Theore'ra 3*7*5 • Suppose 'R 2 holds* r > ^ and that the free 
boundary S is bound'ed away fitsm zero that is there is a positive 
constant ^ such that S ^ v (Bor all € C'0#^l» Then the 
following estimate holds for a backward difference in time^ 


'(3W7*23-) sup {il 'IU —tLl II * “ 0(AT t h D, 


•n 


where the norm iHUlHI ^ is understood in the usual sense 

and 

•= XO# min '(S® 'SCt^)^)) <• 

J^roof* Since sup '|.s^—'W^!| — OCAt+h ) and s > there exists 

'» 
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an s > 0 such that > V-,e =s > 0 (say), for n = 

From ( 3 . 3 . 7 ) ^ ^ (W^+^)2^ we have AtJJ > At for n=0,l,,.,,N 
and hence Atj^ > vj At that is AT^ = o(At) . 

Now is’^-S^l < lS(T’^)^S^|+|s"-s(t") 1 < jejl + K 2 |x^-tJJu 


Further^ 


1 lu’^-u?l 1 


h*‘ 1 -nv ilu(x^)..u?n +i iu^-u(x’^)ll 

H^(q'^) ^ h^(q'^) 


(q”) 


< llu -zril^ + K^iX -tJJi 


Thus* 


(3*7*24) sup 1 +lS*^~s”j} < K(K2,K22HAt4-h^'^^~^’ 

n Ii (q ) 


+ sup !x^—x?l}» 


Now we estimate sup IX^-X^I • Prom (3»2*12) for t = t^ 


(3^7*25) <3^. x’^ =t dj. x’^ - + (s’^'^^) 


n+l % 2 


Subtracting (3*3*7) from (3»7*25^# we get 


(3.7*26) dt(x’^~t“) ^ « l^j + eJ^^(s^'^W‘^^)* 

't 


Since* 




and X t V"-^j > < K '^'= ■'■'tt* ' 2,.-a .n+X> 

n=sO ^ Ln4-1 naO L (t *t ) 

tc 


< K(K«),1JshJJ <5 

- 2 V. 
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we have^Summing on n ^ using -the estimate (3*7*1S) 

and. multiplying ty ^t^ the following estimate 

(3*7* 27) < K(K2)(At + h^) . 

Taking supnenrum over all J/ ^ J 5. n)^ substituting this in 

(3*7*24)# we got the desired result as At, = 0(At)# 

n 

Gotpllg-by 3*7 *6# Suppose that all the assumptions in the 
, o 

Theorem 3*7*5 hold and that Sj.^ satisfies inverse property* 
Then*the following estimate for At = o(h) 

( 3*7*28) sup { I 1 o ^»r, J = 0(h~^ AT. + h^^) 

n ^ ) 

holds* 

Proof * Now, 

< lie’ll!, + HU^-U(T“)n „ 

< lle’^llg + K(K2) 

From the Corollary 3*7*4, aTj^ = 0(At) and the estimate { 3*7*27), 
we get the desired estimate (3*7*28) because ^ o(h)* 



CHAPTER 4 


STEFAN PROBLEM WITH A QUASILINEAR PARABOLIC EQUATION 

IN NON DIVERGENCE FORM 

4«1 Introduction » 

In this Chapter, we deal with a single phase Stefan 
problem with a quasilinear parabolic equation in non-divergence 
form* As usual, we use a coordinate transformation to fix the 
domain and then apply H^-Galerkln procedure* Optimal error 
estimates in L as well as H -norms are obtained for this 
problem* 

In Section 2, the discription of the problem and 

1 

transformed problem are given# The weak formulation and H — 
Galerkin procedure are discussed in Section 3* Section 4 deals 
with an auxiliary projection and some ^proximation Lernnas* 
Optimal error estimates in L^H^ and H^-norms for continuous 
time Galerkin ^proximations are derived, assuming existence 
of the approximate solution in Section 5* Finally in Section 6 
the question of global existence and uniqueness of the Galerkin 
approximation is discussed. 

4*2 P acobl em Description and. Domain Fixing • 

We now consider the following nonlinear Stefan problem 
in a variable domain oC*) x (0,T^] as giv® in Ch^ter 3. 

Find a pair u . UCy.V. s = S(t) such that g satisfies 
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(4.2#l) U^-a{U)Uyy = 0, (y,T) e q(t) x (0,T^] 
with initial and boundary conditions 
(4*2«2) u(y,o) = g(y)# y 8 i; 

(4«2«3) T>0, 

u(s(t),t) = 0, 

and S(t), the free boiindary satisfies 

(4*2#4) S^- = -Uy(s(’T).T)^ T > 0 

with S(0) =: 1* 

Here, the following conditions are assumed on a( •), g and on the 
pair £U, S} and are termed collectively as Condition 

CONDITION B, 

( i) For p 8 R, a(p) > a, where a is a positive constant* 

q 

(ii) For p 8 R, a(p) e C (r) and there is a common bound > 0 
such that 

lal, Ic^pl and 

( iii) The initial function g is sufficiently smooth and satisfies 
the conpatibility condition that is 

gy(o) = g(l) - 0« 

(iv) The problem (4#2#l) - (4*2.4) has a unique solution- 

Further, asstime that the solution {U,S} of (4.2#l) - (4*2.4) 
satisfies the following regularity condition : 
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ubl“’(o,Tq/h^'^^(q(t))) n l“(o,t^;h2(q(t))) 

n n w^'"(o,Tq/w^'®"(q(t))) 

l'*(0,Tq/ W^'* yo-(T))), and 
S e ^‘•'““(O^Tq)# 

Let K 2 be the bound for the functions in above mentioned norms* 

In order to fix the free boundary, we use below a co- 
ordinate transformation 

( 4 * 2 - 5 ) X = s"^(T)y. 

Further, a time scale transformation 

-0 

(4.2.6) t = t('c) = ; s dT' 

o 

is used to decouple the resulting transformed system. 
Following 2*3, an routine calculation shows that the function 
u(x,t) = U(y,i) satisfies 

(4.2.7) Ut"^^^^^xx “ ^ ® I, t e (0 ,t] ; 

(4.2.8) u(x,0) « g(x), x e 1/ 

(4.2*9) ^^(Oj-b) = u(l,t) = 0 , t > 0 

and the function s(t) = sCc) satisfies 

(4.2.10) ^ = -u^(l)s, t > 0 


with 


s(0) = 1. 
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Here t = T corresponds to t= T^, Note that all the regularity 
assxomptions for £U,s} carry over to £u, s} with the bound say, 

¥-2 and the new regularity assumptions be collectively called 
• Purther, the integral (4 •2*6) can be rewritten as 

(4 •2*11 ) ^ = s^(t) with t(o) k o» 

i 

4*3 The Weak Formulation and H ^Galerkin Procedure • 

Consider the i^ace : 

H^d) = {V e H^(I) : V Co) = v(l) * 0} • 

«irv 

The wealc solution of (4*2*7) - (4*2*9) is given by 

o 2 

(4*3*1) ® ^ ' 

H ~Gal erhin procedure * Let c H ( I) be a finite dimensional 

O-r O O 

subspace belonging to regular family* Let satisfy the 

approximation property (l*3*6) with k = 2 and inverse property 

Vi O 1 

(1*3*9)* Now we call u : (0 ,t] Sj^ an H -Galerkin approximation 
of u if it satisfies 

(4*3*2) (^tx'^x^'*'^^^^^^^xx'^xx^ = ^x^^^^^'^xx^' ^ ® ^h 

and the initial condition 

(4*3*3) u^(x,0) = Qj^ g(x), 

o 

where is an appropriate projection of u onto at t = 0 
to be defined later* 

Further, the Galerkin approximations Sj^ and of s and 
respectively are given by 
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(4*3*4) 

and 

(4.3.5) 


dt 




®h 


with s^(0) 


1 



with Tj^(O) » 0* 


4 .4 Some Approximation Leirmas . 

Set 

( 4 . 4 . 1 ) A(uyv,w) = ^ ® 

2 

V and w e H • 

The boundedness and Garding type inequality for A can be 
established (see^ Lernna 2*6.1) by standard arguments* 

Lemma 4*4.1. For u e , v and w e H^(l)/ 


( 4 * 4 . 2 ) lA(uyv/w)l 

and 

( 4 * 4 * 3 ) A(u/v,v) > a 1 Iv^^l 1^ - Pllv^l 1 

where a and P are constants* but M* P may d^end on I lu^I 1^^* 
Define 

(4.4.4) Ap(u/v,w) = A(u;v,w) + 

o 2 

Note that Ap(u|*,*) is coercive in H ,that is 

<10 2 

( 4 . 4 . 5 ) Ap(u;v,v) > a Hv^Il * 

Let u e be an approximation of u with respect to 
the form An t 

to O 

( 4 . 4 . 6 ) Ap(u/u~u,X) a 0, X e 



Now, an application of La:^Mllgram Theoran. (Thaoraa 1.3.3) ahows 
the existence and uniqueness of u, defined In ( 4 . 4 . 6 ). 

ConsMer 


(4-4#7) L*(u)<P (a(u){p )+u (l)(x<p ) - 

9x ^ ^ 3CX X 

For ^ e lHi), define <P e by 


VL e 8 ^( 1 ). 


( 4 •4<»8 ) L*' (u)(p = ^ X e I ^ 




= 0 * 


Then for v e H^(l) we get 


(4*4#9) (v/L''(u)<?>) s= Ap(u/v,<P). 

Thus# defining D(l*) as 

D(L'^) = £<p e n §2 : ^^( 1 ) = 

we have from the positivity and boundedness of that at least 
a weak solution <P e D(L*) of (4«4*8)# for each f e exists* 
Further, by the regularity Theorem (Theorem 1*3*4), we get 


(4*4*10) H<PI I 4 < CqI 

where depends on u and its derivatives* 

Let V = u-^* We now need to obtain some estimates of 
V and its ten^joral derivatives for our future use* The 
following Leiima proves very convenient for our purpose* 

Leittna 4*4. 2* Let § e H^( l) and satisfy 
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In order to get an L —estimate, we follow here Aubin— Nitsche' s 

duality arguments* For f e L^(l), define (p e D(L*) by (4*4*8) • 

Oo 

Multiply both the sides of (4*4*8) by § to obtain for u e H , 


(4*4*16) (f/^f) = (f,L‘''"(u)<P) = Ap(u;§,<P) 

= [An(u;$,<P-X) + F(X-<P) + F((P)] 


< Ml I $ J I inf 1 I<P-X1 I „+M. inf 1 1 <P-X1 1 „+M 1 1 <Pl 1 

XX Q « JL Q A ^ 

XeSj^ Xes^ 

2 


4 


< [Ml I iK^h^ + M^K^h^ + M2 ] 11 <PI 14* 


From (4*4*10), (4*4*14) and (4*4*16), we obtain the required 
estimate (4*4*15)* 

The next Lemma contains the error analysis of 7) and 

Lemma 4*4.3* For t 6 [o,t], the following estimates hold for 

T) and 7)^ 

(4*4.17) ^ ^4 ^ 

and 

(4.4.18) <K 5 h”-''(IIUtll:„+"’^'U.J= 0 ,l ,2 and ?<m<«l, 

Here K, and K. o« poaltlve constants depadlng on parameter 

^ ^ . -v 

e5(pressed through the following e35)ressions #4 4 o'" 3 

and Kg = K5(K^,K^,K3,K^,IIUtll^2,~ 

Proof. Put « = V and P B 0 in the previous Lemma 4.4.2 to get 

I ■ V" - ’^3 ^ 

XtSh 
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< inf 1 1 u-v 1 1 2 # for v = X + u e s. 

< W""2||ull„. 

For Tie IITIII < llfi^ii and llTi^ll < llTlj^ll and hence the 

result (4 *4 #17), for j = 2# 

Similarly, we get the following estimate for 1 1 HI I , 

1 |r)l 1 < K.h"'nul 1 ^# 

“•4 m 

Consequently, the estionate for llnilj_ follows from the interpola- 
tion inequality, 

I ini < I iT)i 1 inf 

In order to estimate T)^, we differentiate (4 #4 #6) with respect 
to t and obtain 

(4.4.19) Ap(u/VX) = 

Identifying the right hand side of (4*4*19) with P(X), we see 
that for <p e H^( l) 

1P(<P)I < Kg iitp^fi, 

where Kg deponds on and 1 Itig 1 

Purther, for (P e D(L*) we get on integration by parts 

(4.4.20) P((P) = )-u^^(l)(Tl,(x(p^)^) 
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02 

For u s H and = 0^ the second term in the right hand 

XssO ^ 

side of (4 *4 *20) vanishes* So* 

whcuTQ K, = 

Thus* Lemma 4*4*2 is applicable with M, = K^HT? !] and M sK-linil 
and we get the desired estimate (4«.4<*18) for j «= 0,2* For j = 
as usual we make use of the interpolation inequality to get the 
estjunate (4%4*18)* 

We Shall also need later the following estimate for 

Lemma 4*4*4ii There is a constant Kg » Kg^a.,K^*iy!ljlK^) such that 
for 2 < m < r+1// 

( 4 V 4 •* 21 ) m ^ 1 ) II < Kgh^^ i| .1 uH :i 

Proofv Define an auxiliary fuasi'Ction <P e iD as a solution o^E 


L*(u)<P » 0, X'e Jjf 


(p ! 

XXXi 


X=sO 


= 0 


■<P I: 

XX, 


= U 


x*l 


Multiplying by ^ the first equation and integrating by parts, 
we obtain 

.A^^xufV/fP-^h, >X e \ 
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< M 11T)|| inf I l<p-xlU 

* O 2 

xssh 

< 1 lul l<Pl 1^. 

Hence the result follows# 

^ *5 A Priori Error Estimates for Continuous Time Galerk in 
Appro ximat ion ♦ - . . 

Throughout this Section, we assume that there are 

positive constants K* and h^ such that a Galerkin approximation 
h ^ 

u 6 in (4 *3 *2) exists and satisfies 

(4*5*l) I lu^I I 00 <5 ^ 0 < h < h-, 

L (H^) “ " ° 

where u“(x/0) is defined as Qj^g, satisfying 

(4.5.2) Ap(g/g~Qj^g,X) = 0, for X 6 S^. 

Clearly, u^(x/0) s u(x,o). 

Let I = u^-u and e = u-u^ =17 

c*3 h. 

Theorem 4.5 #1# Suppose n = u-u satisfies (4.4.6) and u , the 
Galerkin approximation of u is defined by (4.3.2) with given 
as in ( 4 . 5 . 2 ). Further, assume that (4.5 •!) holds# Then there 
is a constant Kg = Kg(a,P,K* K 3 ^,K 4 ,K 5 and Kg) such that for m>4, 

( 4 . 5 . 3 ) 


Proof . Front (4.4.6) and (4.3.l) with v =X, we get 
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Subtracting this from (4.3.2), we obtain 

( 4 .5 .4 ) (S X)-Ap(u;u,X )=(h^^X^)~p(7i^,X^)+p(g 

But, 

(4.5.5) Ap(uV^X)-Ap(u,u,X)=(a(u*')^^X^)+([a(uh).a(^)]~^ 

' 1 >^>+ V 1 ^ ' »4Xxx>-5x‘ 1 > ‘ ’^'Xxx'+ <=< 5 x^>. 

Prom ( 4 . 5 . 4 )— (4 *5 *5 ) with X =s ^ we get on Integrating by parts 
respect to x the two terms on the right hand side of 

(4.5 -4) 

2 ie' '^x"^+<®<'^''>2xx'2xx>=-‘VSxx>+<<’’'^xxWl'<='5x'^xx> 

+([a(u)-a(u‘')] 

Using a( .) > a, (4.5.1) and r^lacing u by u— 7], we obtain 

(4.5.6) ||gllg^ll^+allSxx'‘^<“^t*' llg^II+PJIhll 

+K^\\t^\i^ +k2(ii 77J1 co+ngiJ «> ) 

L L 

+K^K2( 1 1 7)1 1+ 1 is I I ) I I l+K*lh^( 1)1 I \t^\ 1 
+K*lg^(l)l 

1/ 2 o 

Since 12 (l)l < I 12^1 IS^^l 1 ^or t e H^, applying Young^s 

SL^ S 2 

inequality for the lest term and the inequality a3D< b / 
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a^b > Oj e > 0 for the remaining terms in (4.5*6) to get using 

I l<PI I eo < 1 |(P„I 1, for (P e §2 
L ^ 

+ 1 17)1 I ^+l'n^(l)| 2)+K(K2,e)l |T)1 I jl 17)1 I 2+K(K2/e) 1 |T|l f 2 

I ISj,! I VK(K^,K2,#,e) lU^ll2. 


Now with appropriate choice of 6, K^Q(e) can be made less than 
or equal to oc/2. With this choice of 6| we get by integrating 
with respect to 't' and using Gronwall's Lemma (Leirma 1*3*9) 
the following^ 

I lS^Il^(t)+a ; I 1^ dt' < K(Kj^,K2*K*/p);(||n^li2+llT,|l2 

o o 

+ It)^(l)l2+| l7)l l2| iTjl |2) at'* 

Prom (4.4*17), (4.4*18) and (4*4*2l) with 2(n>-2) > m and 
2m~3 > m that is m > 4, we get the desired estimate (4*5*3). 

Corollary 4*5*2. Let all the assumptions of the Theorem 4-5*1 

o 

hold and the finite dimensional ^ace satisfy the inverse 
property ( 1*3*9). Then there exists a constant ■^1 depending 
on such that for r+1 > m > 4, 


(4.5.7) nsii 
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Pi 2 °S. the estimate (4^.3) ana IISI, < nj^ll for J 4 § 

we get, ^ 


By inverse property (1*3.9), we have 


L^(H™)' 


‘ *^xx’ ' oo o < h"*^! IS 1 I 

L (l2) - ° =< ^12) 

Hence the result (4^5.7) • 

Prom Theorem 4.5.1, Corollary 4.5.2, Lemaa 4.4.3 and triangle 
inequality we get the following Theorem. 


4.5.3. Let the solution u of (4.2.7) - (4.2*9) satisfy 

the regularity hypothesis Rj_. Further, suppose that there are 

positive constants h^ and , (i^ > 2K^) such that an approxiinete 
In o 

solution u e in (4.3.2) satisfying (4.5 *1) exists In 
I X (0,Tj for 0 < h < h^. Then the following estimates hold 
for r > 3/ 


( 4 #5 .8 ) ^ ^ ® ^ ^ ®0/ i \ ^ 2 ^ ^ j “ *^#1/ 2, 

L (H^) 

where ^^2 ” ^12^^2^^4 

Besides, for sufficiently small h and r > 3 
(4.5.9) i lu^l I 00 o < 2K«< 

and conseguently, can be chosen Independent of K*. 

Proof . The estimates (4.5*8) for j = 0,1,2 are immediate from 
the Theorem 4*5 .1^ Corollary 4.5.2 and Lemma 4.4.3 by triangle 
inequality. 
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To prove (4 *5 #9), we note 


I fu^l 1 


2 ^ - ' 0 + lleil „ n 

L (H^) l“(h2) l’”(h2) 


< Ka + K^a 


r-l 


< 2Ka/ for sufficiently small h and r > 3* 


We are now looking for approximations of U and S 
where the pair {U,S} is the solution of (4#2*l) - (4 •2*4)* The 
Galerkin approximatiorsU^ and Sj^ are given by 

(4.5.10) U^(y, t) = u^{x,t); 


(4.5.11) = Sh^t), 


where 


y = Sj^(t)x*' 

(4.5-12) 

T 


and are given by (4.3.4), (4-3*5) respectively. 


Theorem 4-5.4. Suppose that the condition B and the regularity 
hypothesis are satisfied. Then the following estimates hold 

for r > 3. 


(4-5.13) 

(4-5.14) 


I IS-S^l 1 

1 I 


l“’(o,t^) 

Ao,t^) 


oCh"-^) . 

0(h”'^) 


and 
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(4.5.15) inu-u^iii a 

L (OiT^;H^(oCc))) 

where 1 I I .1 I I is defined as in (2.7.23) - (2*7.24). 
Proof* From (4.2*10) and (4*344)) we have 


\s-s^\<S (\Vyil)\ + \l^l)\) Isl dt' +;iu^(l)i is-sj^l dt' 


t 
S < 

o o 

An application of Gronwall'g Lemma (Leirnia 1.3.9) and the 
estimates (4.4*2l)/ (4.5*3) form a r+1 give 


(4.5-16) 1 I K(K3)t 1 


< K(Kj) 

< for r > 3. 

The estimate (4.5-13) is iirmediate from (4.5*16), if we note 
that 

Further, the estimate (4*5.14) follows from (4.2.11), (4-3.5) 
and ( 4 .5 *16) . 

Finally/ since 

we obtain the required estimate (4.5.15). 

Restart . For optimal rate of convergence in i^-norm. we can 
use finite elements oonsisttag cl-splines of degree r > 3. 
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4*6 Global Existence and Uniqueness of the Ga l erkin 
AppjTOxijTici'bjLon » ' ^ ■ " 

Now we consider the problem of existence of the 

Vi 

Galerkin approximation u in the domain of existence of u* 
Towards this end, let us recall (4*5 *4) and note 

Ap(u*’;u'’,)0-Ap(uj'u,x) = Ap(u,i:,x)+([a(u^)-a(u)] 

Firom the above we get/ 


(4*6,1) + %Cu;g/X) = 

+ ([a(u)-a(u’^)] ^xx^Xxx^-V^^^^x'^xx^ 

+ 2 j 5 .(i)(xu^#X^jj), 

But, 

(4*6,2) a(u)-a(u^) » a^e « e d| # 

Replacing u^ by u-e in (4*6*l) with (4*6*2)/ we have 

(4.6.3) 


o 

‘ - ■n^(l)(x(u^-e3,),Xxx)+£x'^^'=^’"x-®x>'Xxx' 

Substitute e by E(x,t), where E e H . Then we get, 

(4.6.4) ^ = -(nt,X3tx)+(>t''.X^>+^^x'’^x’ 
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1 

-(/ 

o 


a.(U"«^E)/_ y \ jt. / 


+ S^(l)(x(u^-.E^),X^^) i 

which is a linear ordinary differential equation in t • Therefore, 
for any E = E(x/t) there exists a unique solution t of (4 *6 #4) 
with 


(4*6*5) S (x,0) = 0 

in tho interval (0,T)* 

The equation ( 4 * 6 * 4 ) defines an operator^ such that 5 s=^(e), 

O 2 

for each E e h • Since e =sT) -2, therefore 
(4*6'*6) e » T) — ^j?(E), for E E 

To show the existence of a solution in (4»3#2)/ we need 
to show that tho operator equation (4*6*€) has a fixed point* 

In othorwords, wo are looking for an e(E) such that 


e(E) 


E* 


Theorem 4*6*1* Suppose that the finite element space 

satisfies inverse property (l»3*'9) ^ i® “the unique solution 

of (4.2*7) .* ( 4 . 2 . 9 )* Further, let the regularity condition 

be satisfied. Then for some 6 > Q# there exists a solution 

"h o ih 

U e Sj^ .of (4.3*’2) satisfying ililu-u T/H^'(l)) ~ * 

Proof . Set X. * t '( 4 *'6*4) to get 
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+ ivi)i (iiUxii+iiViiiUjjxii+diUjjii+iiE^n) 

IJx'l’l "Sxx"' 

Using I5^(l)i < applying Young's inequaility 

2 

for the last term and the inequality ab < -^ + | b^,a,b>0,e>o 
for the remaining terms vre get# 

I \2al I ^‘< K^^(s)l I \K(K^,p,K2/e)( I 1 ^+1 I’ll 1 

L 

(1+1 lEl l2)+l'»lxCl)l^(l+ilEll2))+K(PrK^/K2,e)(l+|lEl l|) 

CO . 

Choosing e appropriately so that 2oc = K^^(e), integrating 
with rospoct to t and there after applying Gronwall's Lemma 
(Loirma 1*3*9) we get 

1 I2l 1 J(t)< K(K^,K 2 ,P )e:>p [K(P,K^,K 2 )t(l+MEl *^ 00 ^ ^^2) { ll\l 

/ 

+ (l’l^(l)l^+ ll7llli)(l+llEll2 )}. 

^ L (H'^) 

Prom Lemmas 4*4*3 ” 4*4*4#it follows that 

(4.6.7) 

where = Kj^jClC^.Kj.K^.Kj.Kg.P and 


Thus, we have 
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C4.6-8) 4 IISII 


L (H^) 


l“(h2) 


< 1 1^1 1 


- 1 , 


For I 1 El I 


CO o + h *1 1 gl I - 

L (h2) O ^ l“(h1) 

~ ^ froin (4 #6 #7) — (4#6^8)|- we get 
1 I el 1 „ < K 

L~(h2) - 16 


where K^g depends on and 6 

Therefore, for sufficiently small h 


1 1 ell ^ < 6 * 

L (H^) - 

Now, an application of Schauder' s fixed point theorem (Theorem 
l*3«ll) guarantees the existaice of an E such that e s= E, which 
is a solution of the operator equation (4#6#6)# 

The uniqueness can be proved following the arguments 
of Chapter 2 » We formalize the above in the form of a Theorem# 

Theorem 4#6#2# Let u satisfy the regularity hypothesis and 

h ® 

let K > 0# Then there exists one and only one solution u e 

of (4.3 #2) in the ball {llu-u^ll « < K} , for sufficiently 

L (H^) 


small h and r > 3* 



CHAPTER 5 


SINGLE PHASE SEMILINEAR STEPAN PROBLIM 


5 Int roduct ion • 

In this chapter we consider Galerkin methods for 
approximating the solution of a single phase semilinear Stefan 
problem for continuous time^ backward difference in time as well 
as extrapolated Crank Nicolson schemes# Optimal rates of 
convergence in and H^-norms are derived for these 

cases* 

In Section 2, the main problem, the basic assumptions 
and the cmestion of straightening the free boundary are discussed. 
The weak formulation, Galerkln procedures and some approximation 
L^as are described In Section 3. Section 4 deals with optimal 
error estimates In and HUon® for the continuous 

time Galerkln spproxlmatlon. A prion error estimates for t 
discrete time Galerkln procedures are derived in Section S. 

5,2 Problem ynrnulat l on and Basle Assumptio n s , 

A single phase semUlnear Stefan problem Is stat 
follows = .ind a pair . = U(v.r), S a S(T,1 of functions 

such that G satisfies the aemlllnear parabolic eguatrcn. 


ti _ U = f(u), e QX ' 

yy 


C 5 » 2irl ) 

subject to the Initial and boundary conditions. 
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(5#2-2) U(y^O) = g(y), Y £ 1/ 

Uy(0,T) = 0 , 

(5*2#3) T > 0/ 

U(s(t),T) = 0, 

where the free boundary SCt) is governed by 


(5*2*4) + Uy(S(T),T) =: 0, T > 0 

with initial condition S(o) = 1* 

We assume that f( *) is a bounded smooth fiinction on R, 

and that a pair {U, S} of sufficiently smooth solutions exists 

and is unicjue* 

We organize our regularity hypotheses on the solution 
{U^S} of (5»2#l) " (5*2*4) according to the results in vAiich 
they are used* We assume 


R^# USl'^CO^Tq/H^^CqCi: ))) n L““(0,T^/W^'“’(0(Tr))) 

n w^<2(o,Tq;h^^(q('c)) ) n W-^^CO/T^/h^CqCt))) 
s e w^'''"(o,Tq); 

R 2 . ueRj^ n w^'^Co^Tq/w'-'^CqCt^))) n w^' ^'( o,T^;L^(Q('r))), 


s e ^ n w^'^(o,T^); 


rv> 

U 8 R 2 


nw2'2(o,vwi'“(!i(t))) n w2'2(o,i^;h2(o(t))) 

nW^'^CO^T /L^( q{ '*^)) ) / 


s e R 2 n w^'^(o,Tq)* 
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Let be a boTOd of the fttotlons in all the nor™ of the 
spaces in and R^, 

We now introduce new co-ordinates x and t, connected 
with y and f by the relation (2-3.1) and (2.3.2), If we 
write U(y,T) = u(x,t) and S(t) = s(t) In the transfonned 
variable^ then the problem (5*2«l) - (5*2*4) is transformed into 
the following coupled system of equations : 

(5*2*5) ■Ujj.(l)Ujj.+s^(t)f(u), X e I, t e (o,t]* 

(5*2*6) u(x#o) s g(35)/ X e 1/ 

(5*2*7) u^(0,t) = u(l/t) = 0, t > 0 

and 

(5*2*8) || = -u^(l^t)s, t > 0 

with s(o) = 1* 

The relationship of t with t can be rewritten as (2*3*12)* 

Note that the regularity hypotheses ^^*^2 ^3 solution 

£U, S} can be carried out to the solution £u/s} and call these 
R^fR^ and R^ respectively with a bound (say)* 

5*3 Galerkin Procedure and Some Approximation Lemmas * 

Define H^( l) as in 2*4. If u and f(u) are sufficiently 
regular*’ than u satisfies the following weak form of (5*2*5) 

(5-3.1) (nt;^v^)+(u^^v^)=u^Cl)(«i^v^)-s 2 ct)(f(u),y^),veg 2 


with 


u(x,o) « g(x), X e I* 
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Let Sj^ c 8 ^( 1 ) be a finite dlmenadonal subspaoe 
satisfying both approximation property C 1 . 3 .S) for k = 2 and 
the inverse property ( 1 , 3 . 10 ), 

o\\ consider an I^-Galerkin procedure for ^proximating 

the pair {u,s} of (5.2*5) ^ (5.2,8), Denote the ^proximation 

o 

of u by u : (0 ,t] -♦ and the ^proximation of s by s. , where 

h “ 

the pair {u /Sj^} is defined by the relations : 

(5.3.2) ^ 

With u’^(x,0) = Q^g(x), 

where is the appropriate projection onto and 
ds, , 

(5*3,3) = -u^(l)sj^, 

with ® 

h 

Once the pair {u /Sj^} is known the approximation of t is 
given by ( 2 . 4 , 5 ), 


Next, we define a fully discrete Galerkin proced-ure 

based on linearized modification of the backward differencing 

in time. Let At > 0 , N = T/At 8 and t^ = nAt,ns=0,l, 2 , •♦*,N. 

Let (p” = <p(x,t") and d,. = (<p’^‘*'^-<p")/At, Denote the 

o 1 NO 

approximation of u by Z : {t ,t , ,,,,t } ■* snd the approxi- 
mation of s by W. Assuming that and are known, we detei> 
mine and as follows : 


(5.3.4) (at2S,X^)+(^^X^)=^(l)(x2P+^X^)-(«”)^(f(z”),X^), 


XS 
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and 

(5-3^5) d^W^ = 


Both these discretizations are 0(At) in time# A priori error 
estimates will be presented for and s^~W^ in Section 5* 

Here we remarlc that (5 *3 #4) ~ (5*3#5) with given Z° and W° has 
solutions for sufficiently small At# 

For higher accuracy in time, we sh^ll consider the 
following variant of the Crarik-Nicolson Galerkin scheme# Let 

qjn+l/2 _ ( . ^ 1 t Now replace (5»3#4) and (5#3#5) by 


(5.3.6) 

E(w^)’^(f(Ezi^),X^^), X e 

and 

(5*3.7) d^w”^ = 

Where EZ^.S/Z Z^^l/Z Z^^ and 2,, 

The present scheme gives rise to a system of linear algebraic 
equations in z" and (n > 3) and is solvable for sufficiently 
small At, provided Z°/Z^W^ and z^ are given# We note that 
a starting procedure is needed tc define 2^,2^ as well as 
which will retain the overall aooiracy of the method, A 
predictor-oorrector version of (5.3.6) and (5.3-7) «111 anfflce 


for our purpose# 
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The aneiysls proceeds, followlpg Chapter 2, using an 

auxiliary projection. Per this, consider the following bilinear 
form : 

A(u;v,v/) = for u e W^'““/v,w e H^, 

Let P be chosen sufficiently large that the bilinear form 
Ap(u/v,w) = A(u;v,w) + P(v^,w ) 

X X 

Oo 

is coercive on I) that is there is a positive constant a 
such that Ap(u;v,v) > al Iv^l 1^. Let u e be the Galerkin 
approximation to u with respect to the form Apt 

(5 #3 *3) Ap(u;u-U/X) = 0^ X e 

for each t 8 (0 ,t] # The existence and uniqueness of u in 
(5*3*8) can be easily shown using Lax-Milgram Theorem (Theorem 
1*3*3)« Let h « u-u* We now state the following estimates 
for Tj without proof* The proof can be obtained following 2*6 
or 3*4* 

Lertma 5»3*1# There exists a constant such 

that for 2 < m < r+1 and j =0,1,2 

(5*3*9) 1 Ihl I j < K 3 h®^^ 

(5*3*10) 1 Ih^-I 1 j < K3 h® ^ (I lul + I lUt-l 

(5.3-11) lyi)! < "“I'm' 

Further, we need an estimate for for our future use* 
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Following Douglas, Dupont and Wahlbin [l9 ] , we prova tha 
following Lemma* 

Lorctna 5.3.2. Suppose u e L"(o,T;w”'l'T),then there is a 
constant = K^(P,a,K ,K„,| lull _ ) suoh that for 

2 < m < r+1 

(5*3.12) \\\\\ ^ 

L 

P roof * lor u 6 H (l)^ leb u s Sj^be the Ritz— approximation 
defined by 

= 0, X e 

Then, for sufficiently smooth u, 

1 1 (u-*u)j^l I 00 ^ ^4 2 < m ^ r+1. 

L 

^ A ^ 

The difference (p =s e obeys 

Ap(u;<P,X) = -0^(1) (x(u-u)j^/Xxx^ + ' 

The choice X = <P gives 

al 1 <P^^1 I ^ < K( I lu^l 1^<^P) ll(u-u)^ll^„ '^^xx^‘* 

Since, 

00 *^ iKu-u) li 00 + 

^ L L ^ 

and 

1 1 ‘I’xl I » "«’xx' ' i lKu-fi)^l l^„ , 


the inequality (5*3.12) is proved. 
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We note that \inder the hypothesis (5#3*9) and 
( 5 * 3 * 12 ), it is easy to show that there exists constant K_ 
such that 


(5*3*1 3 ) 


lu 


'X* * 00 -5 >). ^ • 

^ L (L ) l2(h2) - 5 


5 *4 A Priori Error Estimates for Continuous Time Galerkin 
App ro ximat io ns * 

In this Section, we first develop a priori error bounds- 
for the error e =: (u-u^) and e^ = (s~s^)* Let t - u^*-u, h =. u-u, 
then e = 'n~S* In order to maintain a unifonn degree of 
approximation, we define Qj^g = u(x, 0 ), where u is the projection 
of u onto defined through (5*3*8)«f Thus S(xjtO) a 0* 

Wc now compare u^ and u, defined in (5*3*2) and (5*3*8) 
respectively* We make the following additional assureption on 

Assume that there exists a positive constant K*(> 2 IC 5 ) such 

that 


(5*4#l) 




Prom this as a Corollary a bound for * * ®h* ’ 
found out, using (5*4*l)* 

Indeed from (5*3*3) we have, 

ls>,(t)l < 1 + / ‘®h' 

^ o 

using Gronwall's Le»ma (Lenma 1.3.9) ana (5.4.1), we get 

IlShH „ < ew (K*I) < K(K*). 

^ L ( 0 ,T) 


(5*4*2) 
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Thoorem 5*4*1. Assume the regularity condition Rj_ and (5 *3*13)7 
(‘^.4.1) to hold. Then there existfe Kg=Kg(a,P,K^^K^/K 2 /K 2 ,K^/K'’‘’) 
such that 


(5*4.3) IlSll oo i+lleill 00 +pll2il o - < K.h^'^^yfor r> 3. 
L (H-^) L (0,T) L^H^)" 

Further, if h is sufficiently small, then Kg can be chosen 

independent of K’’*’. 

Proof * Prom (5*3*2), (5«3*l) and (5*3*8), we have 

(5.4.4) 

( 1 )( ai^.X^)+s^C f(u),X^)-s^( f (u*’) ,X^) . 

Subtracting (5*3*3) from (5*2*8) and multiplying the resulting 
equation by e^, we get 

(5.4.5) 11^101.1^ =-[>^U)e2+Vl)se^-£^C3.)seJ . 

Setting X = S In (5.4.4) and combining with (5.4.5), it follows 
that 

(5.4*6) ^ I > ^+le^l^)+Ap(u;^,S)=:(T)^3^Sx^-P(^x'^x^'*'^‘ ’ ^x‘ ' 

* 

+T1^( 1 ) ( xx'~5 x' ^ ® ^ ‘ ' 

5 )+s2( f (u)-f(n*') , I ) I ® 1 ' V =®1 


+ se^' 


+-0 X the first two terms on 

Integrating by parts with re^ 
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the right hand side of ^ -i . 

uc ot «,5.4.6; and applying Saequality I of 

Chapter 1 to all terms except the fifth and the last terms, we 

obtain 


( 5 *4 .7) I I I ' ®1 1 I 2xx' * ^176 1 1 1 2 +k^(p,k^,K2,k-^8 ) 

( I I\l 1^+1 IVl I Vlr)^(l)l^)+Kg(p,K^,K2,K(K«‘),K*,e) 

(Iis^l|2+le^i2)+|yi)(xuh^g^),+jg^(l) 

To estxrnate tha lQ.st but one teinti on the right hand side of 
( 5 * 4 * 7 )/ we note that 12 (l)l < / 2 I 12 11$ I 1^'^^ 


an application of Yomg's inequality that is ab < 


Now 




1 . X 


q ' 


” + q = 1/ a,b > Of where 6 is appropriately dhosen, gives 

( 5 . 4 * 8 ) Ux^lHxa^/2xx)l <K'^ 12xU)l l^xx" 

< V2 K* I 12^1 1 12^J 


< 4 ~ I I I 23 oc^ I ^/£or p=4/q=4/3. 


We next use 12v(l^^ obtain a bound of the last 


term 


■XX 


K, 


(5.4*9) 1 23^(1) se^l < -| le^^l +®>I?xx’* " 


Combining (5.4*7) - (5*4*9)/ we obtain 

(5.4.10) (1 12^1 l^+le^l^)+al 12^x1 i^<(se+y 4 ^ » 2 xx' * ^ 



+K^(P,K^,Kj,K*,e)( I ID^I |\| ITII I ^+H1 ^(i)|2)+K3Cp,Kj^,K2, 
K(k'''),K*,s,6)(i|2^1|2+ |ej^|2). 

Now choosing e and 6 appropriately so that (se + 3/4 6^^^) is 
less than cx and integrating with respect to ^t* from 0 to t/ 
we get 

t ^ 

(5*4^11) 1 1 S(t) 1 1^+1 ej_(t)l^+a;i 121 j^dt' < 2K„ 

o ‘ ” o 

t iT)^(l)l^) dt^+2Kg ; (11211^+ le^J^) dt' , 

where we used 1I2II < 112x11 ^ * - * * ^xx* ^ ^ ^ 

In view of the estiinates(5*3*9 )jf (5*3*10) for j = 0 and (5*3»ll)/ 
Gronwall's Lentna (Lemma 1*3#9) yields 

I 1 2(t)l l^+lej^(t)l^+a ; 1 12 ll2dt'<2K3K^ e3?)(2K^t){h^^^^ H 

Now, taking supremum over all t e (0,T), we get the required 
estimate (5-4.3) for 2(i>l) > i+l or equivalently, for r > 3. 

To complete our arguments, we must show that Kg can be chosen 
independent of for sufficiently small h- We use the inverse 
hypothesis (1.3.10), (5.4*3) and (5.3.13) and obtain 
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< Kg + K^Kg h‘*V2 


Clearly, for sufficiently small h 


(5*4.12) 


lul 


^ 00 oo 

L (L ) 


< 2Kg < k\ 


Thus, Kg can be chosen independent of 


From Theorem 5*4*1, Lenmas 5*3.1, 5*3,2 and the inverse hypothesis 
(1*3 * 10 ), \ie get the following Theorem* 


5 #4^2# Let all the assunptions of the Theorem 5#4#1 hold 

along with the inverse hypothesis (1*3*10) for Then there 

n 

is a constant = K^Q(K^,K 2 ,K 4 /Kg) such that for r > 3, 

(5.4*13) 00, j, ss j = 0,1, 

L (0,T/HJ(I)) 

and 

(5*4*14) lie, 1 I 

^ L (0,T) ~ 

If additionally, u e l'"(0,t;vP*'^'*), then 


(5*4*15) 


Hell 


l“(o,t;w^''*(i)) “ * 


Proof. Since IlSll ^ ^ < M511 i we have by (5.4.3) of 

L (L^) -L (H^) _ 

Theorem 5.4*1, combined with (5*3*9) for j = 0 the estimate 
(5*4*13) for the case j « 0. 

Similarly for j ;= 1, we conclude (5*4*13) from (5*4*3), (5*3*9) 
and the triangle inequality r. 

To coinpleto the proof of the estimate (5»4-15), we Invoke Inverse 
hypothesis Cl. 3.10) for 2 aid the relation (5.3.12). 
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Lastly, (5*4»14) follows dlrert-1.r -p 

airectiy from the estimate (5^4.3). 

To prove an optimal order of oonvargenoe to H^-norm, 
we state the following Theory, 

Thaorem 5 *4 *3# Unde>r- i 

er laentical hypotheses as in Theorem 5 *4 *2, 

there exists a constant ^K,Aa,P,K ,K K K k k #k ^ 
such that 


( 5 *4*16 ) 


®* * 00 , 2 < 

L (0,T/H^(I)) “ 11 


Proof* If choice X = St is made in (5*4^), then 


(5.4.17) 


‘x^^xxf 




* ^ ^2 ■*■ ^3 ■*■ ^4 + ^5’ 

Now for Integrating hy parts with respect to x we get 

\ = '^(l^Sx(l)S^(l)~uJj(l)((S^+xS^),S^)* 

Then using the inverse hypothesis (l*3#10), the estimates (5*4»3) 
and (5»4»l) we obtain 

(5.4*18) ILl < lu^l)l Wl^W oo + 

L 

112^11) 

< K(K^,l<^,e)[h**^llSxJ tVj 
+ 3 eIISxt^'^ 
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< K(K^,Kg,Ki^e )h^^+K(K»,K ,e ) 1 1 g 1 ] 2+ 36 1 1 g | [ 2, 

XX Xt 

For Ig/ we integrate by parts with respect to x, and get 

SO/ 

(5.4.19) jl 2 l < K 5 I >ell^i.,coV’'^^^'*^xfc“ + 2K5l|ell^j,^J 12^1 1 

< K(K ,K , 8 )(h^^llel |2 +l|e )|2 J+ 2 ellS . 11 ^ 

^ YT^ ^ 

< K(KQ,Kg,K^Q,e)h^^^V2) + 2e 

The following bounds for and are easily obtained 

(5.4.20) II3I M- II4I < K(p,K3,e) h^^ + aellg^ll^. 

Finally, I3 can be written as 

= s^(f(u)-f(u'^),SjQ;^) - ej(ej^-2s)(f(u*'>, 

Now integrating by parts takes the form 

I' = -s^((f^^(u)-f^(u’^))uJ^ + f^(u^) V^xt^* 

The bound of is then given by 

(5.4-21) 11^1 <K(K^,K2/e) llell^t 26 N^xtJ}^ 

< K(K^,K2,K^0'®^ 

Similarly/ 
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Tho bound for I” is given sj follows t 

2 


11^1 

< Iej^l^+lej^|2)+|e^j4^jg^j2j j 


+ 46 1 12^1 |2 j 

Using the estimate ( 5 * 4 * 15 ), we have 


(5*4.22) 1I"I < K(K^,K^,K2,K^q,kV) + 4e!lS^ll2. 

Combining (5»4*17) ~ (5*4*22) and choosing e so that the term 
1 I I ^ from both the sides get cancelled, we obtain 

(5.4.23) - "^'V^*'"^ 3 tx"^■*■*'^'’'VW‘S'V'^lo' 

Now, an application of Gronwall^s Lerrma (Letrma 1*3*9), Lemma 
5.3*1 and triangle inequality give the estimate (5*4*16). 

But herb the constant K depends on K^. An argument similar 
to the one in Theorem 5*4.1 shows that Kg is independent of K*', 
for sufficiently small h. 

Now, the Galerhin approximations of U(y,T) and S(t) are defined 
by 

U^(y,'^) « u^(x,t) 

( 5 * 4 * 24 ) , V 1 

and Sj^(T) = Sj^(t) respectively, 

wht’re 

y « S^X 

T^(t) 


( 5 .4 . 25 ) 


T tz 



and are given by (5.3#3), (2*4.5)- 

Theorems 5 *4 -2 and 5-4-3 lead to the following result- 


5»4#4# Suppose that holds and that Sj^ satisfies the 
inverse hypothesis (l.3»10)» Then the following estiniates hold 


( 5 *4 - 26 ) 

1 IS-Sj^l 


= 0(h”’^ 



L (0,T^) 


( 5 -4 • 27 ) 

Iit-V 


= 0(h"l 

-3 

L (0,T^) 



ind 


( 5-4.2B ) 


IU-U“| 1 1 _ » 0,1,2. 


l“(o,t^;h^{'£5(‘^))) 

If in addition, U e L** ( 0 ,T then 


( 5 -4 • 29 ) 


I I lU-U^l i I -I ^ *** 0(h*^) • 

L“’(0,T^/wl'nQC‘t))) 


Here 


iAo,T^;w^'P(q(t))) 


is defined by 


(5-4-30) V''P(q('C))^' 

p K 2 or " and SC''^) is given in (2*7*24). 

5 *5 A Priori Estimates fo r the Discrete._Time Galeiri c to 
Approximations - 

In this Section, we first develop e priori bounds for 
the errors (u"-^) and (s‘’-w”) at discrete time levels t®, where 
Z" and W" are defined by (5.3^) anl (5.3^). ^^e re^ts.obtaln- 

ed will be similar to the corresponding estimates for the 
continuous case in Section S.4. Ma assume here that the 
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regularity hypothesis R 2 of Section 3 holds# Let = s-W, 

^ ~ Z^u and e s u.~ Z =: tj ^ » 

Theorem 5 ♦S • Let Z° be determined such that 

(5*5*1) llZ°-gllj < j = 0,1,2. 

There exist constants and h^ > 

such that if h < h^, At = o(h) and r > 3, then 

( 5 »5 * 2 ) sup( I 1 1 I I ‘ 


5 , 5 . 2 ) sup( 115^1 l^+le)^l^)+p^_^l 1 


Proof* From (5*2*3) and (5*3.5), we have for n = 0 , 1 #. ../N-1, 


,n ds 


(5.5.3) a,e; = a,3''.||| e”«ci)] . 


t.it* 


Next from (5*3.1) and (5.3*8) fort == t^^ we get for n=0, 1, . *,1^1. 


n 


( 5 *5 *4 ) C X):=-( dj.h^, X^)+( d^u^ - gt 


.n+1 


,x^) 




XX 

, rin+l y \ 


subtracting (5*5.4) from (5*3*4), it follows that 


n ‘ X 


(5.5.5) ^ 


n+1 


Xx^ 






-P (71?^ / Xx^'*' ^ ' 
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But^ 


(5 0-6) ) 




XX 


Xj^)+( s''+^-s") ( b^Ks^) ( f ( Z?’) ,X^) 


XX 


+ e”(s”+w”)(f(z”),X^^)* 

Substituting (S*5*6) in (5*^5) with X« and coinbining with 

(5 *5 #3) after multiplying we have for n = 

n ^n+l 


(5.5.7) (atS”,S^^)+<clt4f ^ 


3u 

+(— 2? 
^''at 


n+1 


1 )( 1 


+Ptats^ ( f ( zP ) . J "" )^-<t Aw" ) ( f ( Z” ) ,s 2;' > 




"( l)( 


ntl 


'®1 


Here <(p,W> 

by 


X 1 ^ 

+ <s"+’- 

„„ NOW we aeti^ate the left hand side of (5.5.7> 
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(5.5.3) + A (u''h”-^£ '«•’•) 


- 2At t l^-l I5"l i2)+( le”+'-l2-le”l^)) 


+ a< 1^, 

XX 

ITejd)/ we integrate the first three terms on the right hand side 
of (5»5»7) and obtain 

(5.5.9) 


= -d^u”)+pn“’^^, 


< K(Pfe)( 1 Id^n^l l\p^+l l^)+3el 1 


where 


^n+1 2 a 2 

(5.«^*10) a^*(; ll^(.,tMlldtM2<At X nl^ll^dt' 

at^ 

by Schwarz's inequality. We note that, 

(5.5.U) 


n 9t^ 


(5.5 .11) At S P 
n=0 


n=0 at 




Applying Inequality I of Chapter 1, we see that 
(5.5.12) 
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Similarly, we have 


(5*51.13) gn+1) j 

XX ^ XX 

< K(K^/K 2 ;e)((At)^| l(i^.u^| | \\ |Tl“| \ \\ j 2 ) 

+ 3ell2^1jl2, 

(5*5*14) !(s^+s’^‘^^)At <its^(f(z^),£”+l)|< K(Kj^,K2/e)(At)^ld^s”l2 

and 

(5*5*15) iAh a^uJ(l)(3W^'*'^,S5^'^^)l < K(Kg;.s)(At)^ldtu”(l)l^ 

+ e lis"t^l iK 

XX 

FurtVier, tVie following estimates are routinely deduced applying 
Inequality 3 of Chapter 1 


Lt *5 * 


lu ) l?l”( 1 )( 1< K(S ) 1 1 1 1 2f77“(l ) 1 ^+e J 12 

X X XX ** XX 


n+l j I 2 

XX ' 


/ r. 1,1 n <2s{ „n+l. , 7 n+lf,\ -,n+l Jitl. , 

(5*5*17) Kd^s ~ ^ X ®1 '®1 


n+l 


< (iZ^+^(l)l+l)le”’^^l^ + 


where 


t"! .2. 


^n+1 , 


(5*5.13) °l n ~ 

^n dt t“ ^ 


and 


{ 5 .5 .18 ' ) 1 ) , 1 ' 
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< (1)1 ^+K(K„;e ) I 1 Ve 1 1 1 1 

X A 1 XX 

Lastly, using Ig^(l)| < >/’2 I lg”l we estimate by 

generalized Young's inequality (Inequality III of Chapter l) 

L 



Multiplying each of the estimates (5*5*8) - (5*5*19) by At and 
summing up n = 0^1, 2, * * *, J-l, we obtain 

(5.5*20) i{(llg^ll^-.llg°l l2)+( J^l^-le°l^)}-Ki.At 

< At(l2E+e/2)’^S^ J T® ^ 

n=0 XX 

+K(p;e)At S (a„+a^ _)+K(p,K.,K^/6)At E ((in II 
n=o ^ ns=0 

+ 1 1 d^Tl" I I I I I 1 ^ool’l 1 ) I ) I ^)+K(K^ , Kj, 

jj 

K_/e ) (At ) ^ '^S^( 1 1 1 1 Vi i \s^\ I V I dtu^( 1)1^) 

^ n=0 

+K(K. ,K,,P;e)At E (1+1 IZ^"^^! 1^00^ HsV I ^ 

+ pAtl Igjl I Vk(k.,k /e)At E {(n^l+l)V(i+ 
lzVl)l)He?l V/t K(K /e)(lzj(l)l+l)le^l^'* 

Assume further that 


f 
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(5.5.21) llz’^11 


< K^. 




Wg shall use an induction argument to treat lw”l. As an 
induction hypothesis, assume for sufficiently small h 


( 5 . 5 . 22 ) 


IW^I < 2K2 , 


fo3r n = 1 / 2 • # ♦ * 


Using ( 5 . 5 . 21 ) and (5*5.22) in (5*5.20) , we see that if for 
sufficiently small At the coefficients of le^l^ and 

•ia IKl 

can be made less Idian 1/2 and similarly for an appropriately 


chosen e, the first two terms in the right hand side can be made 

less than cxAt S 1^, then we have 

n=0 ^ 


( 5 . 5 . 23 ) 1 


l^+Je^l^+ ^At S 1 12 ” 1 1 2<K(p)(At)^( 1 1-^1 1 2 
^ ^ n=0 " at^‘ L 


2(j^2) 


+ ll“il )+K(P,K.,K ,K /K*)(h^^^^^ 

dt L^(0,T) 1 ^ J 

+h2(2r^l))+K(K,,K /K )(At)2(l j|H,l l2 
1 2 5 3t 

J*'! 

+ )+K(p^K,,K,/K*)At S ( 1 Ig^l |2+|ejl2)* 

L (0,t) ^ ^ n=0 ^ ^ 

An application of discrete-Gronwall^s Lertma (Leirma 1.3.10) gives 
us the estimate ( 5 . 5 . 2 )/ which shows that the induction hypoth- 
esis ( 5 . 5 . 22 ) holds for n = J/ for sufficioitly small h. However, 
the constant involved in the estimate d^ends on K*. By repea- 
ting the arguments similar to that of Section 4, we easily show 
that this constant is independent of K*, for sufficiently small 
h and At = o(h). 
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Finally, an application of triangle inequality together 
with Theorem 5«5»1, Lerrtna 5*3*1 and the fact that I I SI 1 < i IS 11 
and I I < I I, for g 6 yields the following Theorem# 

5*5#2# Let all the assumptions of Theoron 5»5*1 hold* 
Then there exist constants ^o ^ ° 

such that if 0 < h < h^^At * o(h) and r > 3 , 


(5#5*24) sup ( I lei i+hl 1 elf.+ le- I ) < K, .(At + h^"'’^)* 
n 1 1 — 14 

t 

Further, choosing X = At d^S^ in (5 *5 *5) and applying Theorems 
S-S*!/ 5 *5* 2 and Lemma 5 >3*1, one can easily prove the following 
Theorem* 


Theorem 5*5*3» Suppose that holds, r > 3 and that the 
space time discretizations satisfy the relation At = oCh)* 
Then there exists a constant 

such that for h sufficiently small the' following holds : 


( 5 •6*25) sup 1 1 el I - < K + h^ ^) * 


n 


15’ 


We shall now present the error analysis for the extra- 
polated Crank-Nicolson Galerkin scheme* 


Theorem 5#? *4. Suppose that the pair £u, s} satisfies the 
regularity condition R 3 and (5.3#13) holds. Further, suppose 
that Z°, Z^,W^ and Z^ are chosen to satisfy 


(5.5.26) l+le^l+l l^+IK^I I ^+At( I 1+115^=11) 


< Kj^g{(At)^ + . 
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Then there exists constant K. _ssK (paK K K k k tv ) 

ependent of h and At such that for h sufficiently small and 
At 5s o(h) we have 


(5.5.27) sup (||e||+ le^l H-hllel l^)< (it)2 + r > 3. 

*fc 

Proof. Subtracting ( 5 . 3 . 7 ) from ( 5 . 2 . 8 ) fort = 
multiplying the resulting equation by eJ+V2^ 

(5.5.28) 

_2b+l/2(j^)j^+V2 ^+V2)^ 

where 


(5-5.29) 


°l,n = '^t® - cE 


^n+1/2 


From (5.3»6) and (5-3-1)/ (5*3.8) for t = we obtain for 

^ > 2 

(5.5.30) (at^SrX^)+(i;""/^X^)=(<ith^X^)-(4t1^- ^(t'’+V2),^) 

+ [{ E 1 ) ( xZ^+V2^ )^X^"'^X^ l/t^"^^/^) ( xu^( t’^+^Z^) ^ 


+{s^(t’^'^^Z^)(£(u(t’^'^^Z^)),X )-.E(w2)^(f(EZ^),X )}. 

XX ^ XK 

Integrating by parts first/ second and fourth terms on the 
right hand side of (5 *5 * 30 ) for X = g 11 + 1/2 combining with 
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(:3«5.28), we get on summing for n » 

, J-1 

C5.5-31) 2 { (<L n+1/2 n+i/o <? 

n=2 ^ )+<dte^,e“+V2>3+ 5, IIS^V2, ,2^ j. 

n*2 ^ n»2 ^ 

^ - |^( 2)2 n+1/ 2)^p ( ^n+V 2 g n+1/ 2 ) , 

■^^ 2 * ^)+'^i''[ { E, z”( 1 ) ( ^+V 

n=2 o •*• X 

s2;^))]-f i {s^(t"+l>(f(u(t«V2)) jn+V2)_E(„2)n 
na2 ^ XX 


where. 


and 


XX 

J-1 


(f(Ew“),g^V2)j^ S <0 

n =2 n=2 

s(t'«'V2).2n+V2(2)^l/2) n+1/2^ 

X 1 

/xV(a,il,J«V2)^(„^^_^,^^l/2)^p(^n+V2,jn+l/2) 




J-1 


+ + S <a. ■'■ ^3' 

'’2,n = V” - ir 


cr s 

3, n XX XX 
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The first two terms on the left hand side of (5 *5 *31) become 

(5.5*32) S = 1^-IlS^ll^) 

n ■“ 2 

+ ( le^l ^-le^l^)} . 

We now estimate the terms on the right hand side of (5.5.31)* 
For the first four terms 


(5.5*33) I 




~(<^. 1 < 4e I I ^+K( P,6 )'^E^ 

3,n XX - ^^2 n=2 

( i l^+i •«^2jtn* 1 ^) . 


However, the following holds : 


J-1 , ,, 

( 5 . 5 . 34 ) At E ( 1 la. _l 1 >1 la 1 1 '^) 
n=2 

J-l 




In a similar way, we obtain 


^5*5.35) ^‘-n=2 n=2 ^ 


But, 


(5.5.36) At ^ n=2" L2(t",t”+M 


Now we give an estimate of below 
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(5*5*37) 


whero/ 

(5*5*38) 

(5*5*39) 
( 5 *5 .40 ) 

( 5 *5 *41 ) 


2"^^)-(u^( i,t“+V ^)-Eju5J( 1 ) ) ( ^ 2) 

+u^( l,t"+V2, ( ^(■^n+V2_^( x,t"+V2, ], jW-V 2) 

+ |c ( u;+l( 1 )-u^( 2) ( 2)^.( ^ ) 

-u^( l,t" V 2) ) ( n+V 2) j ^ j ) , ^^n+l/ 2^ ^n+V 2, 

<7E J I I^+K(K2,k /e) z''{ ( 17)"(i)| 2+|I7"'2 (i)|' 

n=2 =► n=2 

+ lla^^^ll^} + q + IJ , 


r-i 

1 


j~i 


B -£ lu^(l,t^+V2)^EV(i)|2 

=2 n=2 ® ^ 

< K(At)3 

hho.Tf^n, 




^ l"(l”)- 


o,n XX — xt * 

E^l Ia» _ 1 1 (^n+l/2) , j 2 

=0 n=2 ^ ^ 


ia2 


< K(At)^ 1 

l2(h^) 
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C5*5*42) 

and 


J—l 
I S 

n=2 ^ 




(5.5.43) 3^' - I S E zS(1)(x2”+V2,j.wV2,|. 

ris— ^ 


Further, we see that using l 9 ^(l)| < \/’2 • ^ 

for (f>e P can be bounded by 


< K(K )^2\ 2i$"-V2(i ) ,+3,2n-,V2(^ ) , j i jgn+1/2, , 

ns=2 ^ XX 

< ''<•^>32^ 1 i5^-y2, ,1/2 

,ljJ^y2||V2),|^Wl/2,,^ 

Applying generalized Young's inequality (inequality III of 
Chapter l), we get 

(5*5*44) < e I I ^+K(Ke,e) E^( 1 lg^~V2| I 2^j jgn-.3/2| ] 2j 

^ ns=2 ^ ^ 



+ellS^/2, j2^g^2!l^V2, j2^j^(j^5*e) ^ ( H 1 ^+] 1 2^'’^ 1 1' 

2 

+ 1 1 ^)<2e E^l 1 ^+s( 1 1 1 ^+I 1 h 

X ^ o XX XX XX 

n=2 

or**l 

+ K(ic/e) E lisjjll^* 

^ n=0 ^ 

Using Inequality I of Chapter 1, we have the following 

(5.5.45) lf< K(e)VlE z”(l)l 2 llS^^/^l l^s'^E^ 

n=2 ® n=^2 
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Further the follovd.ng also holds : 


S^(t’^'^^/2)f(u(t”^l/2))^E(w2)^f( EZ^)=(s2(t”^^/2)„E(s2)”) 


-f ( EZ’^) )+( f (uCt’^‘‘*^/2).f ( Eu’^ ) ) j . 

Hence by the inecguality I of Chapter 1, we have 
(5.5»46) ll^ll K(Kj_;e) S^C l<7g^^l2+Iag^^l2)+K(K^,K2/s)^S^( 

+ IIE£"ll\llo,„ „ll^) + 56 E 

1^/’^ n=2 


where, 

( 5 ♦5 *47 ) 


t 

i 2 ='^ 2 ^ 1 E(w 2 )’^-.E<s 2 )’^ 
n=:2 n=2 


2 < ^E^IECW^+s’^) l2lEe”|2, 

" n=2 


(5 *5 *48 ) 


j—l 


H la 
n-2 



2 




and 


( 5*5 *49 ) 
Further^ 


"'i' n 0 , „ 1 1 1 1 Eu''-u(t"+^ 2) , I 2 <k( A t) 3n n ^ • 

1 °'" n =2 ’ 


^(t"+V 2 ). 2 "V 2 ( 1 )w"+V 2 ^(^^( 1 > > 

.(tntl/2)^.zn+l/2(2) 3n+V2^2^1/2( 5(tn+V2)_3n+V2, 

X 
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therefore using 1 £^(i) 1 < i 1 5^1 1 we get 

( 5 . 5 . 50 ) Jl^l < 2^E^l2^+i/2(i)|2ign+l/2,2^/~^j,gn+l/2, ,2 

^=2 ^ n=2 ^ 


where* 


-t S (lT)"+V2Q)j2^jg(^n+l/^gn+l/2j (l)!^) 

n=2 ^ llgTi 

+ K(K„*e) '^S^le?"‘‘^/2] 2^ 


( 5 . 5 * 51 ) 


J-1 

n=2 


J-1 


(l)l^=E Iu’^‘’'^/^(l)-u (l*t’^'*'^/^)l 2 
n=2 ^ 


X 


< K(At)^i 1 


u. 


tt 


I 1 


2 

L^(W^'“) 


# 


Note that 


(5*5.52) E ls(t^'^^/2)-s^‘*'^/^l ^ < K(At)^ I 1 s. . H^’ 

n=:2 “ ^ Lno*T) 

Combining (5*5*3l)-(5.5.52) and multiplying the resulting 

equation by At* we get 


(5.5.53) t( 1 IS’^I I ^)+At E^l ! ^<20e At E^l i ^ 

n=2 n=:Z 

o'*** ^ 

+ K(P*K.*K.,K-;e )za; E ( ! Id^n^l 1 ^+1 I ^ 

^ ^ ^ n=2 

+ I jETpl l^+|7?^'^^/^(l)l^+l7)JJ(l)l\ln”~^(l)l^ 


X X X 

j 

+ m^’”^(l)J^)+K(K^*K2*Kg;6) At S^(At)^ H rgl 1 ^ 


J-l 


n, , 2 


+ K(P,K^*K2*K5?e)(At)^+K(K5/e) At ^ 
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+K(e) AtlEg^"^(l)l2i |gJ| j2^j^ 

^ n=2 ^ ^ 

+K At(l+l^"V2(^) j2) j Jj2_|,j^ At,'^2^(l+lZ^‘‘'^/^(l)l^) lejl^ 

n=2 ^ 

+K(K ) A/s^(lEW^l2+i)(|eyi2+,en~l,2)^g 

n=2 i 1 XX 



Now assume further that for any K'^ > 2K^ 

(5.5*54) IW^“^1/I IZ^I 1 „ < K% for n = 1,2,...^J* 

^ L ~ 

From Lenma 5.3*1 and the estimates (5*5 .53)# (5.5*26) and 
( 5 * 5 * 54 ) it follows that for sufficiently small At, the 
coefficients of and lell^ can be made less than 1/2 and 

IXb «!• 

similarly by choosing e =s we have 

i lS^ll\leJl2<K(p,K^,K2,K3,K^g/K'^){h2^^^W^^^^^+h2^(At)^ 

+(At)^}+K(Kj_,K2/K^/K*) At E^( 11 2^1 1 ^+1 e^l ^) ♦ 

Since At = o(h) and r > 3/ an application of discrete Gronwall's 
Lemma (Leirtna 1*3*10) gives 


(5-5*55) li^^li’^+le'^l^<K (P/K^r^KgfK^g, 


K*){ ( At)'^+h^^ h . 


Following the technique used in Theorem 5*4*1, it is easy to 
show that the constant K involved in (5-5*55) is Independent 
of for sufficiently small h. Further, using 1 U H < ! I 
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a.or S e H , Lemma 5*3»1 and triangle inequality, we obtain the 
desired result* 

SjSSSudS* prove (5*5*27), we need (S*5»26)* For this purpose, 
a predictorv-corrector version of (5*3.6) and (5.3*7) will suffice 
as a starting procedure of sufficient accuracy* Since the 
techniques are similar to those presented in Theorem 5*5*1 
and are similar to those presented in Douglas et al [l7] * 

modified appropriately for H^-Galerkin procedure we shall not 
present them here* Further choosing X = At d^S^ in (5 *5 #30) 
applying Lemma 5.3*1 one gets estimate of e in H norms at each 
discrete time level t^. 


We are now looking for a discrete-time Galerkin 
approximation of {U,S}, where {U, S} is the solution of the 
original problem (5*2.1) - (5.2.4). The approximation 
of {U,S} is given hv 


( 5 *5 *S7 ) 3 .nd 

where y^ = x W^ and is given by either of the two relations 
(5*5*58) d^T^ = (W^'*'^)^ 

or ■ ' 

(5.5*59) ^ "^h " 

Here corresponds to t’^ under the approximate transfermation* 



148 


Consequently, the time mesh in corresponding to a uniform 
mesh At in t will not in general be uniform* 

Let ATj^ s t? and AX^ = max Ai:?* 

0<n<N 

Finally, we establish the following Theorem for the discrete- 
time error estimates U^-uJ and where s U(t^), s’^ss(t^) 

and T corresponds to t^ in the original transformation# 


5»5#5* Suppose that R 2 holds, r > 3 and that the free 
boundairy S is bounded away from zero that is there is a positive 
constant v such that S > v > 0 for all t e [o,T], Then the 
following estimate holds for a backward difference in time 
Galerkin approximation 


(5#5-60) sup { 1 1 +is^-s5ju = 0(At + h^^”^),j=0,l, 2, 

0<n<N ^ ^ 


where the norm 11*11 is imderstood in the usual sense and 

h^Cq’") 

q” = ( 0 ,min (S^,s^,S(t^)) ) . 

Proof# Since sup Is^— = 0 (At + and ^ there 

n 

is an e > 0 such that ^ > v-e = > 0 (say), for n=!0,l, 2, ** •,N« 

From (5#5#58) ^ we have At|J > At, for 

n = 0,1, and hence AT^ > At that is ATj^ = 0(At)# 

Now applying the arguments in Theorem 3*7«5 and Theorem 5*5»2 - 
5 # 2# 3 we get the estimates (5 *5 # 60 )# 


For an e3q)olated Crank-Nicolson Galerkin scheme, we have the 
following results* 
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Theorem 5.5*6* Suppose that is valid along with the 
assumptions in Theorem 5 - 5 . 5 . Then 

( 5 .5 *61 ) sup { 1 1 I . ^ + 1 } = 0£ (At ) , j=0, 1 

n h-J(q") ^ 

holds* 

Proof* It is easy to see that >0 and A^^^ = 0 (At), 

Now 

< \e^\ + 

and 

liu^-u"ll . ^ < lie’ll!. + K j = 0 , 1 . 

H^(Q ) j ^ n 

Thus^ 

(5.5.62) sup { llu’^-uJJnj+lS^-S^n < K(K 2 ,K^ 7 ){(At)\h^^'"^’ 

+ sup 
n 

To estimate we have from (2.3*12) for t = 

and ( 5 . 5 * 59 ) 

^ h ^'^n+l /2 . 

S ince# 

lat - •'h' i A ‘ - ■'S' 

)l<K(K 2 )(At) 

n=l ^ i^n+1/2 

and < «Kj) At, 

n=l 
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we have summing on n = using the estimate (5*5*27) 

and multiplying by At, the following estimates 

(5*5.63) < K(K2) {(At)\h^'^^} + 

Assuming (5*5*26), = 0{(At)^ + and taking supreitum 

over all J (0 < J < n) in (5*5*63) we get 

(5*5*64) 1 K(K2) {(At)^ + * 

Now combining (5*5*64) and (5*5*62), we get the required results* 



CHAPTER 6 


extensions and remarks 
6 *1 Int roduct ion » 

In this concluding chapter, we make some informal 
observations regarding extensions of our analysis to more 
general problems such as the ablation problems and fluid 
filtration in porous media etc*, and discuss various modifica- 
tions which might improve upon our analysis in certain directions* 
More overawe present a critical assessment of the problems and 
the methods, used for solving them* Further, we present a 
short discussion on superconvergence phenomena# 

6*2 Refinements and Generalizations # 

Throughout the thesis, we deal with the free boundary 
problems either with homogeneous Neumann or, with Dirichlet 
boundary conditions* But the above analysis can be easily 
extended to the free boundary problems even with time dependent 
boundary conditions* For exairple, in Ch^ter 2 if the flux at 
y = 0 is given by 

(6#2.1) Uy(0,T:) = <P(t) , T > 0 

then the corresponding condition (2.3*10) for u can be rewritten 
as 

(6*2*2) = <y(T(t))s, ■ 
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In this way the nonlinear parabolic equation (2*3*3) - (2*3#9) 
and (6*2r2) with u(l) =0 is coupled with two ordinary differen- 
tial equations ( 2*3*11 ) and (2*i3*12)* The error analysis follows 
the same lines as the corresponding time independent case with 
qualitatively the same bounds* Similar analysis is also valid 
for the free boiandary problems with time dependent Dirichlet 
boundary conditions that is 

U(0,T) = (p(T), 

Let Sj^ denote the approximation of the free botindary/ 
then the error s-Sj^ is of the same order as u^(l)— u^Cl)* Since 
X = 1 is always a knot, a superconvergence occurs* Following 
quasiprojection technique for Galerkin method for heat 
equation [ 15 ] , we can achieve a higher order of convergence at 
the knots than the global error for nonlinear Stefan problems* 

Next< a more general single phase Stefan problem s 

= 1^ (a(y,T) ||) + bCy^T) || + c(y,T)u + fCy^t)^ 

with either Neumann or Dirichlet boundary conditions, can be 
tackled using the approach of Chapter 5* 

6*3 Ablation Problems » 

In an ablation problem a solid is melted by a moving 
fluid* It is generally assxmed that the diffusion coefficient 
is so large that = 0 in the liquid. Thus, the heat input 
qCi) from fluid to the solid is prescribed on the unknown melting 
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ary y _ s(t), giving rise to a one phase Stefan problem [38 ]• 
Problem P. Given T > 0 and g(y) for y e I with g^Co) * g(i) = 0 . 

Find £U(y,T),s(T)j such that S(T) > V > 0 for 0 < T < T, a; being 
a constant, 


(6.3*1) 

( 6 . 3 * 2 ) 


(6.3*3) 


( 6 • 3 *4 ) 
and 


S(0) = 1 , 

Ut-- (a(y)Uy)y = 0, T e (o,T] and y e q(t), where 
Q(T) {y . 0 < y < S(t)}^ 

^(y^o) = g(y) for y e I, 

Uy(0,T) =: >f(T) 

for 0 < r < T, 

u(s(t),t:) =0 


(6*3*5) 3 t " ^yj + q(‘’^) forO<T<T. 

y=s(T) 

Here, we make the following assuirptlons on q,a,g and the solution 
{U, S} and collectively call then/ condition B' • 

(6-3*6) CONDITION B* 

( i) The function a( •) is twice continuously differentiable and 
has a uniform boxind such that 


> max (lal, ISyl, la^yl)* 

( i±) There is a positive constant a such that a( •) > a. 

(iii) q(T ) is a bounded smooth function of r such that lq{,iq^I<Kj^ 
( iv) The pair {U, S} is determined uniquely and is sufficiently 


smooth « 
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or simplicity^ we take fCt) = o, although the case fC"^) > 0 
can be treated similarly* 

We state ejjplicitly the regularity hypothesis on the 
solution {U,S} of (6. 3*1) - (6*3*5) by as given in 5 #2. 

a Problem P can easily be transformed into a Problem P 
with a fixed domain by a Landau-»type transformation [383 
y = SCt) X 

(6.3.7) 

-c = t. 


If we write U(y,T) = u(x,t) and S(t:) = s(t), thai by the routine 
calculations of 2»3 we get the following transformed Problem P. 

P roblem P i Find a pair {u, s} satisfying 


(6.3.8) s\i^-(a(xs)u^)^ = -a(s)u^(l)xu^+sq(t)xu^(x/t)el x (0 ,t], 

( 6 . 3 . 9 ) u(xfO) = g(x) for x s 1/ 

( 6 . 3 * 10 ) u (0,t) = u(l,t) * 0 for 0 < t < T 

and 

(6.3*11) ^ = ~a(s)u^(l)s’*Vq(t) for 0 < t < T 

with 

s( 0 ) = 1 * 

Note that the Problem P Involves a nonlinear parabolic problem 
( 6 . 3 . 8 ) - ( 6 . 3 . 10 ), coupled with an initial value problem for 
ordinary differoritial equation ( 6.3*11 ). Further, under the 
transformation the regularity assunptions on {U,S} are automa- 
tically carried over to £u,s} and are denoted collectively by 
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with a, hoMnd. K 2 ( S 3 y ) « 

fonmlation and Galerkin p rocedure (see Chapter l)# Let 
H - £v e H : v^(o) = v(l) ss 0}* By weak solution of (6*3*8) - 
(6*3*10), we mean a function u = u(x,t) such that for each 
fixed t 8 (0 ,t], u(*,t) e H^(l), u(x,o) = g(x) and u satisfies 
the following equation 

(6*3*12) s (^^v^)+( (a(xs)u^)^v^)=a( s)u^(l)(xu^v^) 

-q(t)s(xu^,v^), V e H^(l), 

For each of the values of h in (0,l], let c: be 

a finite dimensional subspace satisj^ing the approximation 

property (1*3*6) for k = 2 and the inverse property (l*3#9)* 

Vi o 

The Galerkin approximation u“ ; (0,t] - of u is defined by 
(6*3*13) 

with 

ujj(x,0) = Qi^g(x), 

o 

where Qj^ is an appropriate projection of g onto to be defined 
later and Sj^, the Galerkin approximation of s is given by 

(6*3*14) ^ = -a(sj^)u^(l) sJJ^+q(t) with Sj^(O) = i* 

The error analysis follows the procedure adopted in Chapters 3-5 
using an auxiliary projection* Let 
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(6-3.15) A(s,u;v,w) , ((a(xs)v^>^w^).a(s)u^(l)(xv^w^) 
+q(t)s(xv^w^) for u e v and w e 

i or the operator A, the boundedness and Garding type inequality 
similar to Chapter 5 hold with the modification that the 
constants M and P depend on Isl along with other terms* 

Let Ap(s,u;v,w) = A(s,u;v,w) + P(v /W )♦ Then we have 
( 6 i» 3 * 16 ) Ap(s,u/v^v) > CL i Iv 1 1 ^ for u e and v s H^* 

T, - O 

For each t > 0^ we define u b Sj^ as the Galerkin approximation 
of u with respect to the form Ap • 

(6«3*17) Ap(s,u/u-u,X) = 0 for X e §j^. 

Since Ap satisfies both boundedness and coercivity properties# 
an application of Lax-Milgram Theorem (Theorem 1»3*3) gives 
unique solution u of (6 #3*17) : 

Let T) = u-u* Following the analysis of Chapter 2 
or 4# one proves the following Leiima* 

Lemma 6*3*1* There exists a constant = K^CM^a^p/K^^K^ and K^) 
such that for j = 0#l/2 and 2 < m < r+l# 

I ir]] 1^. < Ivl t 

1 JT?^1 I < K 3 h’^^( 1 lu^i Ijjj+I lul 

and 

l\(l)J < K 3 
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Assvime that there exists a constant such that 

^^-•3-18) II ^ < K 

L (l2) - 4 

This follows from the previous Lemma. 

est faiat es * We assume that for positive 

(6-3.19) 

This is not a restriction, since s{=s) > v and we shall 
subsequently see that 


' o> = 0(h^'^^). 

L (0,T) 


Further, we assume that there exists a constant K* such that 


oo, 9 < K*, where K* > 2K„. 

L (H^) ~ 2 


(6,3.20) llu^ll 

] 

We shall see later that for small h, this Indeed holds good* 
Now we define the projection Qj^ as follows : 

(6.3»21) Aq( s,u;u— Q^u,X) » 0 for t = 0# 


Clearly, u ^x,0) « u(x,o). 


We use below the standard notations = s— Sj^, e 
t = u^~u and e = 


h 

u«u , 


Theorem 6 .3 .2. Suppose that the regularity condition and 
(6 *3 *18) - (6 .3.20) hold* Further suppose that the projection 
is defined by (6.3*21). Then there is a constant depending 
on M, P,a,i',KQ,Kj_,K 2 ,K^ and K* such that for r > 3, 
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(6#3*22) IlSll 


i,"(h 1) .+nsii 


< K,h 


L-(0,I) - l2(h2) - “S' 

E£22f* From ( 6 . 3 * 12 ), (6»3*13) and (6*3»17), we have 


iH-1 


+ ^ s^-s^) (u^^X^)+( ( a( xs)-a( xsj^) ®“®h^ 

(a^( sx)u^,X^)+s( (a^( xs)-a^( 

-a(sj^)r)x(l)(xu^X^^)+(s-Sj^)q(t)(xuJj,X^3^)-( ( ■ 

-( a( s)-a( sj,)) u^( 1 ) (xu^X^)+a( 1 )( xu^X^) . 
Subtracting (6*3*14) from ( 6 » 3 .ll) and multiplying by e^g we get 

(6.3*24) i |t*eil^=a(sj^)e^(r)^(l)-.g^(l))s"^+a(sj^)u^(l)(ss^)''^ e^ 


+ (a(s) - a(sj^)) ■u^(l) s"^ e^* 

Setting X = 2 in (6»3.23), adding the resulting equation to 
(6«3*24) and following the analysis of Chapter 5 with the 
additional information that s^ > we obtain 

1 iE-(s^(t)ll£xl'^+leift+alUxxll2 < Kg(h2^ 1 1 2+ je^ I 2), 

where Kg depends, in particular upon K*. Then applying Gronwall's 
Lemma (Lemma 1»3»9), the required result follows* 

Corollary 6*3.3. Let all the assunptions of Theorem 6 •3*2 hold. 
Then for r > 3/ 



159 


(6 #3*25) 
Further, if 

n 


satisfies the inverse property (1*3.9) then 


^t.3.26J II^j, < j^r 

L (h2) - ^ 

The following result gives the estimates for e and e^. Since 
e analysis follows exactly the pattern of proof in Chapter 5, 
we state only the Theorem without proof. 

^eor^ 6,3.4. There is a constant depending on and 
such that for r > 3 


(6.3.27) 

and 




(6.3.28) lle.Il < 

^ l'”(o,t) “ ^ 

Besides, for sufficiently small h and r > 3, 

11“*’“ CO, 2^ < »4<K* 

L (H'^) 

and consequently, K,^ can be chosen ind^endent of K% 

Finally, the Galerkin approximations of {U,Sj are defined as 
in (2.7.3) with y = Sj^(t)x and t = where Sj^ and are 
given by (6,3.14) and (2,4*5) respectively. 

Theorem 6,3*5. Suppose that the condition holds, then the 
following estimates are valid for r _> 3, 

(6.3.29) IIS-S^II „ = oCh”^), 
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(6-3*30) 


where j 1 1 I 1 1 


'-3 = 0 , 1 , 2 , 

given by (2.7*23) - { 2 * 1 * 24 :) • 


l”(h^(q(t))) (2.7*23) - { 2 * 1 * 24 ) • 

^ For optimality, the finite dimensional spaces used 
3-J^© ^ "Splines of degree r > 3. 

(ii) It is to be noted that a more general problem like 
the one in Vasilev [62]< where the basic equation is of the 
form 

Ot - ®(U)Uyy = f 

can be analysed by the method of the Chapters 3,5 coupled with 
the above technique* 


6 »4 Fluid Filtration in Porous Media* 


In this Section, we consider a semidiscrete finite 
element approximation to a free boundary problem for the pressure 
head of a compressible fluid flowing through a homogeneous 
porous medium* Ey straightening the free boundary by a 
co-ordinate transformation, a continuous-time Galerkin procedure 
is discussed and a priori estimates are derived* Earlier 
Nitsche [51] initiated the study of error analysis for a linear 
problem, posed by Magenes [41 , Problan III]/ but the present 
one is an e 3 (tension to the nonlinear version of it. 

The mathematical mod^ of the phenomenon can be formally 
stated as follows [ 6 ] : 
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Problem P. 

sC^) > 0 , 0 

Given > 0 , find {U(y,T),s(T)} such that 

< ■" < ^o' 

(6.4.1) 

S(o) = 1 , 

( 6 * 4 * 2) 

U^-(a(U)Uy)y = 0 for (y#T) e q(t) x (0#T^] / 

(6*4.3) 

U(y,0) = g(y) for y e If 

( 6*4*4) 

Uy(O^T) a -.f(T:) for T > 0/ 

(6*4*5) 

and 

U(SCC)^T) a s(T) for T > 0; 

( 6 * 4 * 6 ) 

S_ + a(U)u 1 a 0 for > 0* 


Here f('C) > o and is a smooth function* For siirplicity^ let 
f(T) = 0, but for f(T) > 0, all our analysis can be carried out 
easily with appropriate modifications* 

Assume the same condition and the same regularity hypothesis 
R 2 as in 2*2 for a(*), g and the solution {U,S}. of (6*4*l) - 
( 6 ♦4 ^ 6 ) • 

By a* co***ordinate transforinatlon (2^3^l) a.nd (2»3 j^ 2)/ 
the Problem P is reduced to 

Problem P* If U(y,'c) = u(x.t) and S(t) = s(t), then an easy 
calculation as in 2*3 shows that the pair £u,s} satisfies 

(6*4*7) u^.-(a(u)uj^)^ = -a{u(l))u^(l)xuj^ mix (0 ,t] t 
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u(x#0) = g(x), X e !• 

Uj^(0,t) = 0, t > 0/ 

u(l^t) = s(t)^ t > 0 

+ a(u(l)) u^(l)s = 0^ t > 0 
s(0) » 1, 

where T corresponds to of t • Moreover^ t and are related 
ty ( 2*3«12)* 

Remark * The transformed Problem P involves a nonlinear parabolic 
eqnation (6*4*7) - (6*4*10) in a fixed domain, coupled with two 
ordinary differential equations (6*4 •ll) and (2*3*12)* Further, 
the regularity hypothesis R 2 is easily carried over to the 
solution {u, s} and these regularity condition be termed 'condition 
R 2 ^ with a bound (see 2*3) • 

Weak solution and Galerkin approximation * Let H ( l) s= £v e H (l)t 
v^(o) = 0} * Prom the Problem P, we get u^(l) = -•a{ s ) su^( 1 ) . 

The weak solution u of (6*4*7) — (6 *4*10 ) is given by 

(6..4.12) (a(u)u^)^,v^)+sa(s)u^(l)v^(l) 

= a( s )u^( 1 ) ( xu^ v^) * V B 8^. 

Let Sj^c be a finite dimensional space satisfying 

the approximation property (1*3*6), k = 2 and the inverse 
property (1*3*9)* Now the Galerkin approximation £u of 


(6*4.8) 

( 6 .4 *9 ) 

(6*4*10) 

and 

(6*4*11 ) 
with 
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{u,s} is defied as follows ; Find such that u^:(0,t3 

satisfies 


( 6 » 4 * 13 ) 



= a(s^)u^(l)(x.,h X e 

with u (x/O) = where Qj^ is an appropriate projection 

later and Sj^ satisfies 


( 6 *4*14) 


L(Sh)u“(l)sh with Sj^(o) = 1, 


Further the approximation to t is stated as in (2*4*5)* 

The error analysis proceeds^ following an auxiliary 
projection as defined in Section 3* Here in the present case, 
Ap(s,u/v,w) is given by 


(6 *4 *15) Ap(s^u/V/w) = SL(u)v^)^w^)-a( s)u^(l)(xv^w^) 

+ P(v /W )/ u 6 v,w e H^» 

•At 

Following Chapter 3/ an auxiliary projection as well as the 
related estimates for can be easily given* The estimates S 
in H and H -norms can be deduced using the analysis similar to 
the previous Chapters* But the only term, which requires 
special treatment is s sCsj^) 

.< and IS^(1)I < « and 

using Young's inequality (Inequality II of Chapter l) for p =4, 
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<5 4/3 one gets error estimates in 2 ♦ Now the basic estimates 

of e and e^^ read as follows ; 

6«4*1» Let Ujr s and g be sufficiently snooth and 
be defined as (6»3»2l) for t = 0# Further, let satisfy the 
inverse property (1.3.9). Then the following estimates 


i le! 1 


’ChI) 


+ hi lell 


l“(h2) 


+ 1 1 




(0,T) 


hold for r > 2. 


O(h^) 


Finally, the Galerkin approximations and of U and S 
respectively are defined similar to (2.7.28) and the main resiilt 
concerning the errors, in U-U^, S-S^ and t can be stated as 
a theorem. 


Theorem 6 .4 .2. Under the assuirption of the previous theorem 
and the regularity condition and the following estimates 
hold for r > 2, 


s-s^l 


+ |1T-X^II = OCh'^) 

(0,1^) " h (0,T^) 


and 


11IU-.U^1I1 ^ i - j = 1,2, 

L (0,Tq/H^(Q('c))) 


where 111.11! 


l®"(o,Tq/h^(q(t))) 


is defined as in ( 2.7. 23)“ ( 2.7.24) . 


Remark . Similar results can be derived for Dirichlet boundary 
conditions following the procedxire of Chapter 3 and 5* 
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^•5 Remarks and an Overview * 

In Chapters 2 and 3# optimality in and norms 
only is achieved, and H^-estiraate is established with the help 
of the inverse property (l*3*9)* However, optimal orders of 
convergence in L as well as H'*' and H -norms are derived in 
Chapters 4 and 5* For optimality in the finite dimensional 

O r* 9 

spaces Sj^ belong to a regular ^ family ( for a definition see 

•• 

Chapter l) with at least r > 3* Here we note that H'‘'~Galerkin 
procedure suits well the Stefan problem with a nonlinear 
parabolic equation in non- divergence form, although higher 
regularity is assumed on the solution ( see Chapter 4) • We 
observe in Chapter 5 that unlike its proceding Ch^ters, we 
do not require an Inverse property (l»3#9) and (l*3*10) is good 
enough for our error analysis* 

In the Chapters 2,3 and 4, we have discussed the global 
existence and uniqueness of the Gaierkin approximations and 
consistently used a Landau type transformations in these 
Chapters to fix the domain, where as the time-scale transforma- 
tion helps to decouple the resulting system* This helps to 
solve for u Independent of s* This greatly simplifies our 
analysis* However, the transformed problem as given In Chapter 
5 is a coupled one due to the nonzero source term. Even without 
a time scale transformation,the error estimates can be 
sstabllshea followtog 6.3- This fact is ussful to Sealing with 
fully discrete Gaierkin schemes# 


f 
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We make the followteg observation regarding the 

symmetrization of the operator A in Chapter 2* Note that we are 

really not required to symmetrize since it does not bring any 

2 

extra advantage# Rather in a , Certain sense/it conplicates L 
analysis for the auxiliary projection (see the difficulties 
regarding the choice of the domain of the formal adjoint 
accosicated with syntnetrized operator Ap)# The error analysis 
related to the auxiliary piojection in Chapter 3 is different 
from that of Chapter 2 and is more general in its approach* 

Another interesting observation about Chapter 3 is 
the following : Assuming the existence and uniqueness of the 
Galerkin approximation in a ball 1 lu^l ! «, 2 ^ - K > 0/ 

the error analysis has been derived and later to justify the 
above assuirption the existence (global) and uniqueness has been 
proved in a neighbourhood of the exact solution* We believe, 
this argument can be carried out to study the error analysis 
for a nonltoear parabolic problan by H^-GalarXln procedure- 
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